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Chapter 1 


Functions 


1.1 What is a function? 


Example 1.1. The area A of a circle depends on the radius r of the circle. The rule that 
connects r and A is given by the following equation, 


A=nr’. (1.1) 


With each positive number r, there is associated one value of A and we say that A is a 
function of r, denoted by A = A(r). 


Example 1.2. The human population of the world P depends on the time t. The table 
below gives estimates of the world population at time t, for certain years: 


t (years) 1900 | 1910 | 1950 | 1960 | 1990 2000 | 
P (millions) | 1650 | 1750 | 2650 | 3040 | 5280 6080 | 


For instance, P(1950) = 2,560, 000,000 and P(2000) ~ 6, 080, 000, 000. 


For each value of time t, there is a corresponding value of P, and we say that P is a 
function of t. 


However, in this case we don’t know the rule that connects t and P at the moment. In 
fact, it is one of our tasks to find and understand the rule for a function. 


1.2 The definition of a function 


Here we will look at more precise definition of a function. First let us define a set. 


Definition 1.1. A set is a collection of objects, ‘things’ and states. A ‘piece’ or object 
of a set is called an element of the set. 
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Example 1.3. The following are examples of sets: 


(i) Z = {all whole numbers} = {...,—2,—1,0,1,2,...}, 
(ii) N = {all positive whole numbers} = {1, 2,3,...}, 
(iii) R = {all real numbers}. 


If we have two sets A and B, then a function f : A > B is a rule that sends each element 
x in A to exactly one element called f(x) in B. 


A B 
J 


Figure 1.1: Function f ‘maps’ elements from set A on to elements in set B. 
Many-to-one relationship. 


Here we call the set A the domain of f and B is known as the range of f. Suppose x is 
some element of A which is denoted by x € A, then = is called an independent variable 
and f(a) is called the dependent variable where f(x) € B, i.e. f(a) is an element of B. 


NOTE: f is not a function if it is multi-valued, i.e. if one element from A maps to 
two distinct elements in B. Such a mapping is sometimes called multi-valued function, 
however this is just terminology, strictly by the definition for our purposes f would not 
be a function in this case. 


A B 
f 


Figure 1.2: The above ‘map’ f is not a function, t.e. when it is one-to-many or 
multi-valued. 
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1.3. Representing a function 


Example 1.4. Let S be the function from the real numbers to itself, that is S$: R—R 


given by 
S(a) = 2”. (1.2) 


Here we have defined the function S algebraically using an explicit formula. 


That is, S is the function that squares things (real numbers), i.e. a function performs 
some sort of ‘action’. Here we have described the function verbally! 


We can also represent a function numerically (by table of values): 


Emm ea (ee ea ee 
Co ee ee es 


Finally, a function can also be represented visually, for example by a graph. 


Figure 1.3: Graph displaying the function y = S(x) = 2°. 


If the domain of S' is R, then the range of S$ is R* = {all positive real numbers}. 


Summary: 


A function is a rule, it takes some independent variable (in the domain) for which 
there is some dependent variable (in the range). 


There are four ways to represent function. 


In the following, we take a look at some specific functions.! 


1End Lecture 1. 
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1.3.1 Exponents 


We have seen the function f(x) = x. In general, a function of the form f(x) = x*, where 
a is a constant, is called a power function. For example, 
f(x) = ae a=5, 
i (2) Sa - Gaps ee 


What about f(z) = a"? Can we define a function of this form? Yes, we can, but we need 
to explore the meaning of a”, when x is not a positive integer. 


REVISION: when we wish to multiply a number by itself several times, we make use of 
index or power notation. We have notation for powers: 


x 


a 
a=a-a, a@=a-a-a, a=a-a...a:ad, aE€R, «ceEN. 


Here, a is called the base and x is called the index or power. We also know the following 
properties (laws of exponents) 


deere Sia? so" ek 29 S260 Roa se ae 8 ap. 


DGS Gee | Bo Ske SFB e DR ge a Be ge eae, 


B. ha® ob? (ab) |. e.g, 098? 95982280853 = (2. 8)(2s3)(2+ a) (2 3)3. 


These laws hold for any a € R and z,y EN. 


Now we want to generalise the notation for the power to include whole numbers, fractions 
and irrational numbers, so that a” makes sense for most values of a and all « € R. The 
idea is to make laws of exponents hold generally. 


First, we choose a? = 1 (a # 0) so that 


so that 


aa? =a? =a =1, ora®-a”™ =a" =a) =1, forallneN. 
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al/3 = Ya (the cube root of a, a € R), 


al/" — fq (the nth root of a. If n is even, a > 0; otherwise any a € R is O.K.) 


so that 

1,1 

Pil ee ve een eee Re 
Lipa ct 
Gh velit GiB = gitsts <i) = @, 
n 
——<——— 
1 1 1 n 

ai/®.gi/n.. gli/nm Q@ntat Te Son =o =a 


If x is a rational number, then x = p/q, where p and q are integers and g > 0. Then 
a® = aP/4 = (q)/9)P = (a?)!/4 (order doesn’t matter). 


Therefore, a” makes sense for any rational number 2. 


If x is an irrational number, then we can always find two rational numbers c and d which 
are sufficiently close to x and which satisfy c < x < d. So 


Cen AE EI GOS TE OG: 


It can be shown that there is exactly one number between a°® and a’. We define this 
number as a”. 


Finally, an exponential function can be defined by 
f(z) =a", «ER, 


where a is a positive constant. The domain of f is R and the range is R*. Graph: 


y| 7] 
1 1 
x x 
(a)y=a", a>. (b) y=a", O<a<l. 
Figure 1.4: Comparison of exponent graphs for different values of a. 
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Summary: 


If a and 6 are positive numbers, x and y are any real numbers, then we have 


1 at®t¥ = a®- a’, 


There is also a special exponential function, f = e”, we will investigate this further later 
in the course.” 


?End Lecture 2. 
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1.3.2 Polynomials 


A polynomial is a function P with a general form 


P(x) =a + a2 4 age? +++) + ana”, (1.4) 


where the coefficients a; (i = 0,1,...,) are numbers and n is a non-negative whole num- 
ber. The highest power whose coefficient is not zero is called the degree of the polynomial 
P. 


Example 1.5. 
| P(x) 2 | 3e7+4e4+2] ce | ve | 1-30-4205 | 2t+4 
| Polynomial? | Yes Yes No No Yes Yes 
| Order 0 2 N/A | N/A 3 1 


Importance of polynomials: analytical & computational points of view 


Degree 0: P(x) = ap = aox®, say P(x) = 2. This polynomial is simply a constant. 


yt 
2 
A=R : B= {2} 
+z 
(a) y = P(x) =2. (b) Entire domain mapped to one point. 
Figure 1.5: A straight line parallel to the x-axis. 


Degre 1: P(x) = ax+b,a#0. These are called linear, since the graph of y= ax+bisa 
straight line. The linear equation az + b = 0 has solution « = —6/a. 
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Figure 1.6: Linear graph given by y= ax +b. 


The slope or gradient of y= ax +b is a. It can be worked out as follows: 


change in height 


slope 
4 change in distance 


change in y 
change in x 
(av + b) — (au+ 6) 
v-—U 
= a(v — u) = 78 (1.5) 


VU—U 


Degree 2: P(x) = ax? + be +c, a #0. This is known as a quadratic polynomial. The 
quadratic equation ax? + ba + c = 0 has solutions 
—b+ Vv b2 — 4ac 


2a 


(1.6) 


r= 


Proof. Start by re-arranging the equation and dividing through by a so that 


next we add b?/(2a)? to both sides, hence 


, bP cB 
r+ -x4 = . 
a (2a)? a (2a)? 


Now we can search for a common denominator on the right hand side (RHS) and 
complete the square or factorise on the left hand side (LHS), i.e. 


(» | zm) oH | aaF 
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Taking the square root of both sides we have 


i b 22ac b2 +/—22ac + b2 
f+ =a t = : 
2a (2a)? (2a)? 2a 


and finally re-arranging the above equation (or minus 6/2a from both sides of the 


equation) we get® 
—b+vV b? — 4dac 
2a , 


/ 


Example 1.6. Consider the quadratic polynomial 
P(x) = a7 — 3x + 2. (1.7) 
We can represent P(x) in a different way by factorising it, i.e. 
P(x) = (a — 2)(a@ — 1). (1.8) 


Again, we can represent P(x) in a different way, this time by completing the square. When 
completing the square, we do this based on the value of b as follows. We add and subtract 
(t)* to P(x) such that P(x) = 2? — 32 4 ee 
the equation). Now it is easy to see 2? — 3x + ( 


factorised. Thus we can finally write 


Pays (s- a _ (1.9) 


(2)? + 2 (i.e. we don’t really change 


3 


ay? is the same as (a = aye i.e. it can be 


Now, (1.7), (1.8) and (1.9) are all equivalent and they are each able to give an insight on 
what the graph of the quadratic function P(x) looks like. That is 


(i) Equation (1.7) tells us the graph of P(x) is a “cup” rather than a “cap” since the 
coefficient of x? is positive. Also we can easily see P(0) = 2. 
(ii) Equation (1.8) tells us P(x) =O atx =landx=2 


(iii) Equation (1.9) tells us P(«) is minimal at z = 3 and P (3) =— 


since anything squared is always positive! 


Now that we have some suitable information regarding P(x), we are able to produce an 
informed sketch: 


3 Aside: O = Q.E.D, where Q.E.D = “quod erat demonstratum” which means “which was to be demon- 
strated” in Latin. 
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Figure 1.7: Quadratic graph given by y = x? — 3x 42. 


Degree > 3: Things get a bit more complicated! 


In general, we have the algebraic equation 


ao + az + agz? +--+ + nz = 0, (1.10) 


which has n roots, including real and complex (imaginary numbers, z = a + 13) roots. 


nm=2 we have formulae for roots (quadratics) 
n=3 we have formulae for roots (cubic) 
n=4_ we have formulae for roots (quartics) 


n>4_ No general formulae exist (proven by Evariste Galois) 


But in any case, we may try factorisation to find the roots. If you factorise a polynomial, 
say 
P(a) = (@— 1)(a + 3)(a +4), 
then you can easily solve P(x) = 0, in this case 7] = 1, rg = —3, 23 = —4. 
NOTE: this depends on the property of 0 on the RHS. You can’t easily solve 
(2—1)(e@+3)(@ +4) =1. 
Conversely, if you know that P(a) = 0, then you may factorise P(x) as 
P(x) = (a — a)q(«), 
where q(x) is some polynomial of one degree less than that of P(z). 


Example 1.7. Consider P(x) = «3 — 827 + 19x — 12. We know that x = 1 is a solution 
to P(x) = 0, then it can be shown that 


P(x) = (a — 1)q(z) = (2 — 1)(2? — 7x + 12). 


Here P(x) is a cubic and thus q(x) is a quadratic. 
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Example 1.8. Consider P(x) = x3 — x? — 3x — 1. By observation, we know 


So x, = —1 is a root. Let us write 
P(x) = (2 + 1)(a? + ax +), 


then multiplying the brackets we have 


which should be equivalent to x? — 2? — 3x — 1. Thus, comparing the corresponding 
coefficients we have 


l+a = -1, 
a+b = —38, 
b= =. 


This set of simultaneous equations has the solution 


So we can write 
P(x) = (a +.1)(a? — 22 — 1) 


To find the other two solutions of P(x) = 0, we must set (x? — 22 — 1) = 0 which has 


solution 72,3 = pave = 14+ v2, together with 21 = —1 we have a complete set of solutions 


for P(e) = 0: 


Example 1.9. Consider P(x) = x? + 3a? — 2x — 2. An obvious solution to P(x) = 0 is 
z1 = 1. So we put 


Pt) = (eS 


I 
8 


Comparing corresponding coefficients with our original form of P(x) we gain the following 
simultaneous equations: 


a-1 = 3, 
b-a = -2, 
—b = —-2. 


These have solution b = 2, a = 4 and so we have 


a? + 3x7 — 29 —2 = (a —1)(a? + 4 + 2). 


The solutions to (x? + 4x + 2) = 0 are a23 = —2 + V2, completing the set solutions to 
Bey 0: 


NOTE: In the above examples, the leading coefficient of P(x) i.e. the coefficient of x? is 
equal to one! 


As with many areas of mathematics, there are many ways to tackle a problem. Another 
way to find g(x) given you know some factor of P(x), is called polynomial devision. 


CHAPTER 1. FUNCTIONS 12 


Example 1.10. Consider P(x) = 2° — x? — 32 —1, we know P(—1) = 0. The idea is that 
we “divide” P(x) by the factor (a + 1), like so: 


x? — 27 —1 
z+1) ge —x*-32-1 
—g? —y? 
— 2x7 — 3x 
Qu? + 2x 
—xz-1 
GPL 
0 


Hence, multiplying the quotient by the divisor we have (a+1)(a?—2r—1) = a3—2?—3a-1. 


Summary: 
The highest power of the polynomial is known as the order of the polynomial. 


The roots of a quadratic equation ax? + bx +c, can be found using the formula 


—b+ Vv b2 — 4dac 
2a ; 


t= 


No general formula for n > 4. However roots can be found of higher order polyno- 
mials by factorising first. This can be done by observation, the method of comparing 
coefficients or by long polynomial devision. 


Later we will see that polynomials are functions which are easy to differentiate and 
integrate. 


We can approximate most “nice” functions by polynomials, at least locally i.e. using 
power series.* 


“End Lecture 3. 
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1.3.3. Trigonometric functions 
The measure of an angle 


Degrees: One full turn is 360°. 


Radians: (Usually abbreviated as “rad” or the superscript “° ” meaning “circular mea- 
sure”). Suppose we have a disc of radius 1, we choose a couple of radii and want to 
measure the angle between them. To do so, we measure the length of the arc between the 
radii. Suppose the length is t. This is what we call the size of the angle, i.e. it is another 
measure for the length of the arc. 


A \, 
SO) 


(a) 360° is a full turn. (b) Addition of angles. 


Figure 1.8: Some properties of angles. 


We can see immediately that a full turn is 27rad because a circle of radius 1 has a 
circumference 27. Therefore we have 


lturn = 360° = 27rad, 
1 
5 tun = 180° = mrad. (1.11) 
So, 
1 ° 
lrad = sy 
T 
T 
1° = -=~rad. 1.12 
180" oe 


If the radius is not 1, then you need to take the ratio 


arc length 
——_—— = le (i di i 1.13 

ie angle (in radians) (1.13) 
Radians and degrees have one important property in common: if you follow one angle by 
another then the total angle is the sum e.g. t+ s (see Fig. 1.8(b)). This means that 
addition of angles has a geometric meaning. 
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Trignometric functions: cosine, sine & tangent 


The values cos(#) and sin(@) (usually written cos 6, sin@) are the horizontal and vertical 
coordinates of the point C' (see Fig. 1.9). 


yh 


A i 


O 


cos 0 


Figure 1.9: Geometric definition of cos and sin, circle radius r = OC = 1. 


OA 
cos 9 = oc = OA, (1.14a) 
AC 
ind= ~~=A 1.14b 
sin OG C; ( ) 
AC sind 
tan@ = Oh eae (cos 6 # 0). (1.14c) 


A number of properties are clear: 


1. cos? 6 + sin? 6 =1 


1. 


cos” @ + sin? @ = (46) () _ ACE + BC? _ OC? _ 


OC OC OC? OC? 
Note: AC? + BC? = OC? by Pythagoras’ theorem. 

2. cos and sin are periodic functions with period 27, i.e. for any x, cos(x+27) = cosa, 
sin(z + 27) = sin. 


In general, if f(a +7) = f(x) for all x, then we say that f(x) is periodic function 
with period T. 
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Figure 1.10: Graph of cos @. 
Figure 1.11: Graph of sin @. 
Note: cos : R > [—1, 1] and sin: R > [-1,]]. 
3. cos is even, i.e. cos(—0) = cos. sin is odd, i.e. sin(—0) = — sin. 
In general, if f(—x) = f(x) for all x, then we say f is even. If f(—x) = —f(x) for 


all x, then we say f is odd. 


4. Shift: cos and sin are the same shape but shifted by 7/2, which means 


cos (5 — 0) 
cos (0 = >| 


5. Addition formulae: 


=sin@, or sin (5 — 0) = cos 0, 


=sin@, or sin (0 — = = cos 6. 


cos(@ + ¢) = cos 0 cos ¢ — sin @ sin d 


sin(0 4 


+ ) = sin@cos¢+ cos Osin d 
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6. Double-angle formulae: 
cos(20) = cos(@+6) 
cos” 6 — sin? 0 
cos” @ — (1 — cos” 6) 
= 2cos?9—1 
= 1-2sin’6. 


sin20 = sin(@+0) 
2 sin 6 cos 6. 


7. Half-angle formulae: Let 20 = a, then 


1 
cos a = 2.cos” (5) —-1l1 =) cos? (5) ean 
2 2 2 
1 wd 
cosa = 1 — 2sin? (5) —= = sin? (5) = — 


Exercise 1.1. Try find cos(@ — ¢) and sin(6@ — @) using properties (3) and (5). 


We define tangent (tan) as follows: 


Figure 1.12: Graph of tan@. 


NOTE: tan: R— R. Also, as 0 > +5(2N — 1), tan@ — oo for N €N, this is the same 
as saying, as cos@ > 0, tan@ — oo. We say the function has vertical asymptotes at these 
points. Also note that tan@ has a periodicity of 7! 


The double angle formula for tan can be calculated using the definition as follows: 


tan(@+¢) = eee — “ 


sin cos ¢ + cos 6 sin @ 
cos 6 cos ¢ — sin @ sin d’ 
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dividing by cos @ cos @ we get 


tan@d + tang 


6 = : 
ae ae 1—tané@tan¢@ 


Also, if we divide property (1) on page 14 through by cos? x, we have 


; sin? x 1 


T 
cos? cos? x’ 


which can be written as 


1+ tan? x = sec? x. 


The secant, cosecant and cotangent functions are defined as 


1 1 
sec x = ——, cosecr = ——, cotr@= 
COS X sin x 


tan x 
a 3 : 
v2 1 
1 2 
t 6 
4 v3 
1 SoH 
(a) 1/4 triangle. (b) 1/3, 7/6 triangle. 


Figure 1.13: Some well known results for particular angles can be derived by the above 
triangles for sin, cos and tan. 


Summary: 
1 radian = 180/7, and 1 degree = 7/180. 
The ratio of arc length to radius gives the angle (in radians) 


cos@ and sin@ have a periodicity of 27, whilst tan@ has a periodicity of 7 (convince 
yourself why!). 


The range of cos@ and sin@ is [—1, 1], whilst the range of tan @ is (—co, ov). 
sin @ and cos@ can be obtained by shifting one another by 7/2. 


Take note of the various angle formulae!® 


5Fnd Lecture 4. 


Chapter 2 


Differentiation 


2.1 Rates of change 


Suppose we drive from London to Birmingham (100 miles). We plot a graph of the distance 
traveled against elapsed time. We want to measure how fast we traveled. The average 
speed of the trip is calculated as follows: 


100 miles 


= h. 
2 hrs onaRe 


distance 


DOO ieee a te , 
(miles) H 


i > 
2hrs time 


Figure 2.1: Graph showing distance travelled against time, from town 1 to town 2. 


However, when travelling on the motorway you do not stick to one speed, sometimes you 
do more than 50 mph, sometimes much less. The reading on your speedometer is your 
instantaneous speed. This corresponds to the gradient (slope) “at a point” of the graph 
given in Fig. 2.1. 


In the following, we want to find the slopes of curves. 


18 
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2.1.1 What is the slope of a curve at a point? 


The slope of a curve at a point is the slope of the tangent line at that point. 


x x 


Figure 2.2: Curve y = f(x) with tangent line at (%, 9). 


A straight line touching the edge of a curve (but not cutting across it) is called the tangent 
line at the touching point on the curve. 


We hope that for mathematical curves, we find slopes “algebraically” . 


Example 2.1. Let us start with the example of the curve y = S(x) = x, 


1 1h L 


Figure 2.3: Curve y = S(x) = x? displaying small increment about x = 1. 


Look at the point (1,1) on the curve. We want to understand how this function behaves 
when wx is close to 1. We therefore choose some « close to 1, ie. © = 1+h where h < 1. 
So we have 

SA +h)=(1+h)? =1+2h+h? = S(1)+2h+h?. 
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This tells us two things: 


1. Near 1, the value of S(1) is also near 1, since if h is small, then 2h + h? is also small. 
This is the concept of continuity, which says: 


f(xo) + f(z) as 2 > 2, 
for a continuous function f. 


2. If you increase x from 1 to 1+ h, then you increase the value of S from S(1) to 
S(1)+2h+ h?, which is an increase of approximately 2h, since if h is small then h? 
is very, very small. (e.g. h = 0.001, h? = 0.00001). 


From (2), we know that as you increase the value of x from 1 to 1+ h, S in creases by 
approximately twice as much as h. Therefore the rate of change at 1, which is defined as 
the derivative of the function $(a) = x7, is S’(1) = 2. 


What about when x = c? As x changes from c to c+ h, S(x) = x? changes from c? 


to (eth)? = ce? + 2ch +h? = S(c) + 2ch + h?. The change is roughly 2ch. That is, x 
has changed by h and x? has changed by 2c times as much. So $’(c) = 2c. We say the 
derivative of S at c is 2c. In general, 


So (a) Sor. 


So, we have the idea that the derivative of a function f at x = c is the coefficient of h in 
the expansion of f(c+ h). 


Example 2.2. Let us consider the function q(x) = x°. At  =c +h we have 
g(c+h) = (c+h) = 03 + 37h + 3ch? + AP. 
Therefore, given what we have learnt, we can say 


d(c) =3c? ie. q'(x) = 32”. 


PROBLEM: Suppose we want to examine the function r(x) = 1/x. How do we expand 
r(c+h) =1/(e+h)? 

Before we tackle this problem, let us turn back to Ex. 2.2. First let us write q(c) = c?. 
Now, to get the change from q(c +h) we have 


q(c +h) — q(c) = 37h + 3ch? + h®. 
We now look at the ratio 


changeing — q(e+h)—q(c) _ 3c?h + 3ch? + h3 


= 3c? + 3ch + h?. 
change in x (ec +h)—c¢ h Pte, 


We now examine what happens as h approaches zero. In this case 


lim {3c? + 3ch + h?} = 3c’. 
h-0 


Geometrically, we calculate the slope of the chord joining points on the curve (e.g. ¢(c) 
to g(c+h)). 
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y' 


q(c + h) [Pye ya e.tpeg eo eee eee ete OY f 


Cc ct+h Ge 


Figure 2.4: Graph showing chord (line) joining the points (c,q(c)) and 
(c+h,q(c+h)) on the curve y = q(2). 


The slope of the chord is 
gle + h) — g(c) 
h ? 
and we watch what happens to the slope of the chord as c+ h gets closer and closer to c 
(i.e. h gets smaller and smaller) in the hope that the slope of the chord will approach the 
slope of the tangent line! 


(2.1) 


Example 2.3. Now let us consider the function r(x) = 1/a. In this case we have 


1 1 


cth c« 


r(c+h)—r(c) = 


Now, let us consider the ratio 


ee © - (5 =) = - (<5) = -ae5) 


and as h — 0, we have 


r'(c) Se Le. r' (x) = =o xr #0. 


NOTE: r(x) = 1/2 is not well defined at 2 = 0 and in this case, nor is its derivative. 


1End Lecture 5. 
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Definition 2.1. In general, we define the derivative of a function f at x as 


f(z) = < (f(v)) = lim f(z +h) — f(x) 


h->0 h (2) 


provided that the limit exists. If the limit exists, we say f is differentiable at x. If we 
simply say f is differentiable, we mean f is differentiable at all values of x. In this case, 
f' (x) ts also a function of x. 


Comments: 


1. We interpret the derivative as the instantaneous rate of change, or geometrically as 
the slope of the tangent line. 


2. Equivalently, 


xox Ty — x 


since if you put h = 2 — 2, thn h>0 — > x, « and f(x) = f(a +h). 


Example 2.4. An example of a function which is not differentiable at a certain point: 


jo) =Wi={ ee 


—x «<0 


Figure 2.5: Graph of y = |a|. 


At x = 0, f(z) is continuous but not differentiable, since through the point (0,0), you can 
draw many, many tangent lines. We can also show 


for h>0O, AC eae =" On, 


for h< QO, OFM = 10) ohne L 1, 
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doesn’t exist! Taking the limit from both sides must give the same answer. 


Example 2.5. What is the derivative of x” for a positive whole number n? Let f(x) = x”, 


then 


df _ 4 f(e+h)— f(z) 
ae h 
= lim Ci is 
h-0 h 
; {x” Me na”—!h ae OEM hep? a hey — 7” 
= lim 
h-0 h 
-—1 
= lim {na n(n ) np eae mn} 
h-0 2 
= ne! 
ie. 
f(z) =2", then f’(x)=na"?. (2.4) 


We could now go on and find derivatives of “all” algebraic functions by definition. But it 
is too time consuming and impractical. 


Important: If f’(a@) = 0, then the tangent to the curve f at x = a is parallel to the z-axis. 


(a) Local maximum. (b) Point of inflection. (c) Local minimum. 


Figure 2.6: Different options for when f'(x) = 0. 


Often, f will have a local minimum or maximum at some x = a. 


2.1.2 Rules for differentiation 


Some simple functions: x%, a”, sinx, cosa. 


Complicated functions can be derived from these simple ones, by addition, multiplication 
and composition. 
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Example 2.6. 


‘ de es 
a+’, ga, xsing, «—a* =a? +4 (-24), =—-sina, cos(x”). 
x 


So it is too time consuming to differentiate each individual function we can think of by 
first principle, i.e. using the definition. 


Instead, we want to build a machine to help us differentiate various functions. The machine 
should contain three bits: 


1. Sum rule, 
2. Product rule, 


3. Chain rule. 


The idea of the machine is to tell us how to differentiate functions which are built from 
simpler pieces as long as we know how to differentiate the smaller pieces. 


BONUS: you can still relate the derivative of the entire function to those of the smaller 
pieces, even if you don’t know what the small pieces are. 


The sum rule: 


If f and g are differentiable, then 


d d 


© (F(a) + 92) = (F(@) + £ (G@)) = Fe) +90) 


Example 2.7. Consider the function f(x) = (a? + e then using the above we have 


(x?) < (x*) = 327 + 42°. 


If you repeatedly apply the sum rule, you have 


d d 


© (file) + fale) +--+ fal2)) = © (AiG) + (l@)) + + 2 al). 
The product rule: 
If f and g are differentiable, then 
© (Fa)ala)) = F(@ala) + Fla)al(a). (2.5) 
Example 2.8. 
< [(2?+1)(@?-1)] = 2a(2? — 1) + (a? +182” 


204 — Qe + 3x4 + 32? 
5a* + 3x? — 2a. 
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Here we have put 
f(x)=a2? +1 => fi(x)=22, 


and 
g(a) =29-1 => g(x) = 32”. 


Example 2.9. Consider the derivative of x°, so 


d 


5) > 4 
ey =e) 
= 4°.¢+2*-1 
bat, 
as expected, since 
d 
ee )Sn 


Here we have put 


and 


This shows that differentiation can be approached in different ways, using what you feel 
most confident with.” 


?End Lecture 6. 
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A special case: if g(a) = c, constant, then g(x) = 0, so 


d d ; 
da (CF (2) = eg (F(#)) = ef (2). 


x 


Therefore the sum rule can be generalised as: If fi, fo,..., fn are differentiable and 
@1,42,...,@n, are constants, then 


[an file) + angola) +++ + an fu(2)] = ar fila) + aa fila) +--+ anfi(2) 


Example 2.10. Consider a polynomial of degree n with constant coefficients, i.e. 


p(x) = ay + aya 4 agz? + agz? +--+ an2”, 
then 
p (x) = a1 + 2agxz + Base hee Nanw” |, 


a polynomial of degree n — 1, with constant coefficients. 


The chain rule: 


The chain rule tells us how to differentiate “compositions” of functions. If we have two 
functions f and g, the composition, denoted by fog (name of new function), is the function 
given by 

fog(x) = f(g(x)), (Dog then f). (2.6) 


Figure 2.7: The composition f og first employs g from A to B, then f from B to C. 


Example 2.11. If f(w) = w?+1 and g(u) = Vu, then 
fog(z) = f(g(x)) = f(Va) = (Va)? +1=24+1, 


9° f(x) = 9(f(2)) = g(e* +1) = V2? +1. 
Here, f: R>R+,g:R* —R* and fog:Rt = Rt, gof:R-R*. 


But actually the function h(x) = x +1 can be defined as h : R > R, so be careful when 
generating a function by composition, take note of the difference between h(x) and f og 
in this case. 
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In general, fog#gof. 


Composition can be generalised further for more functions, for example suppose we have 
three functions f, g and h, then 


fogoh(@) = f(g(hw))). 
If f and g are differentiable, then 


—— (f(9(x))) = f'(g(x)) - g(a). (2.7) 


Example 2.12. Consider the function y(x) = (2?+2a)!°. Here we will choose f(w) = w!® 
and g(x) = x? + 22, (so f’(w) = 10w? and g/(x) = 32? + 2). Then 

d d 

qe YO) = GF) 


= “ (2° + 2a)2") 
= f'(9(z))9'(x) 
= 10(x? + 2x)° - (3a? + 2). 


Essentially, what we have done is to substitute g(a) = 2? + 2x in our function for y(x), to 
make the differentiation easier. 


Example 2.13. Consider the function y(2) = 1/3. We know how to differentiate 1/z. 
So let us choose g(x) = x? and f(w) = 1/w. Therefore we have fo g(x) = y(x). We know 
the derivatives of f and g are f’(w) = —1/w? and g’(x) = 327. So 


= (4) = fo@ye) 


x 
ts 2 
~ wp 
3 
_ a 
i.e. we have d 
ae? 3274 


e.g. 
d,: 
ae — 3657364. 
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If n is a negative integer, then m = —n is a positive integer. So 


aed, eye 

ag) = eo 
*(F(9(2)) 
= f'(g(a))g'(2) 


T 

wy 
5) 

“~~ 
a 


Here we have simply chosen g(x) = 2” and f(w) = 1/w, where g/(x) = ma™—! and 


f'(w) = -1/w?. 


Therefore we now know that 


is true for any whole number n € Z. 


What about when n = 1/2 ie. f(x) = x?. How can we differentiate this? First let us 
think about what we know about «2. 


We know that (x2)? =a! 
Let us consider the function g(y) = y? and take the composition of f and g, that is 
1 
go f(x) = 9(f(2)) = (#2) =a. 


In this case f and g are inverse of one another. What does it mean for g to be the inverse 
of f or f to be the inverse of g? 


Aside (NFE): 


If we take a point x = a in the domain of f say, and it takes the value b = f(a) in 
the range (or image). Then the inverse function takes the image point b and sends it 
back to the point x = a. In other words we return ourselves back to where we started. 
There is a well defined rule that goes from a to b and a well defined rule that takes b 
to a. 


The inverse is usually written as f~', this is just notation. 


f:AvB 


ay = 
f--(f(x)) =a then ernie 
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Figure 2.8: The composition fo f~! = f~!of takes you back to where you 
started! 


Definition 2.2. The function f—! is called the inverse function for a well defined 
function f then 


fof @=f ofa) =z. (2.8) 


Not all functions possess inverses. For example the function f(x) = c, constant. 


Figure 2.9: Multi-valued functions do not have inverses (obvious from picture). 


A function is called a one-to-one function if it never takes the same value twice, that 
is 


f(a1) # f(v2) whenever 2x1 4 22. 


Figure 2.10: One-to-one functions have inverses (obvious from picture). 


Only for a one-to-one function f, then f~! exists. 
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Example 2.14. f(x) = x”, if f : (—o0, 00) > [0,00) then f~! does not exist since for 
each f(a), there are two possible x values corresponding to it. 


Figure 2.11: Graph showing the relationship between y = x? and it’s inverse. 


(End NFE) 


If it was that f : [0,00) — [0,00) ie. considering the positive x-axis only, then f~! 
exists. Which in this case we call f~!(x) = g(x) = \/2, since 


And we also know that 


9(f (x) = g(a?) = V2? = 2, 


i.e. f and g are inverse of one another. 


If you draw a function and its inverse on the same coordinate plane, they must be 
symmetrical about the line y = x. Why? Rotate the xry-plane 90° anticlockwise and 
then flip across the vertical axis. This is because we want the inverse function f~! 
who’s range is the domain of f and vice-versa. Also, equivalent to switch « © y in 
y = f(x), rearranging the equation for y to give the inverse. 


I 
— 
8 
Nir 
Ww 
bo 
| 
8 
or 
=a 
cS 
n 


Now let us return to take the derivative of the function g(f(zx)) 
d d 


de [9(f(z))] = an *) 
I (fla) F(x) = 1. 
Now since g/(y) = 2y, we have 
2- f(a) f(a) =1. 


Finally rearranging for f’(x) we see that 
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i.e. the method is the same for integer n (as shown on pg. 23), in words, “bring the power 


down, reduce the power by one”.? 


3nd Lecture 7. 
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So we have proven that 


a 
dx 


is true when n is any integer and also when n = 1/2. 


=a, 


It turns out that for all possible powers of p: 


d be 
a) = par? 


WARNING: If p is a fraction, such as p = 1/2, then we require x > 0. We don’t want to 
take the root of negative numbers, i.e. \/-m, where m > 0. 


Example 2.15. Consider the functions g(x) = «+ — x”, f(u) = u-3 and so g(«) = 
Ag? — 22, fi(u) = —4ur8. The composition gives f o g(a) = (x* — a?)-3, So 
d _4 
a @*- 28] = F@))s'@) 
= —-(x4 — x?)~3(4e3 — 22) 


Generalisation: Suppose we have three functions f, g and h. Then the composition is 


fogoh(a) = f(g(Aw))), 


and it’s derivative is 


ole 
acd 
o——™ 
S 
Feat 
pe 
——~ 
8 
— 
= 
co 
Fee 
Pan} 
Fact 
8 
— 
= 
— 
8 
—" 


ae (x3 +0)? +1] — 


h(x) = 23 + 2, = hi (x) = 3a? +1, 


gu) =ub +1 = gw) = 50, 
fw)swi = fw) = Gud 


The composition works as 
f(g(h(x))) = f(g(a3 +a) = f((2? +2)2 +1) = (x3 +a)? + 1 oe 


1 
Example 2.17. Now consider the example where y(«) = (2° + x)? + x| *. It is similar 


66,99 


to the above example, however, when replacing “1” by “x” in the square bracket, finding a 
nice composition becomes a little harder. Instead of trying to work out what this function 
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is as a composition, we simply apply the chain rule using the method “differentiate outer 
bracket, work inwards”. When the derivative is performed on y(x) it follows the steps: 


‘ (3 +2) +2]! = ; («3 a)? z| st < (a a)? x| 
= ; (x3 x)? z\ AG ay Sa’ + x) 4 1 
= ; («3 x)? | AG x)~2 (322 + 2) + 1 


Note that in addition to the chain rule, the sum rule has also been applied. 


Finally, there is another way of writing the chain rule (or other derivatives). Let w = g(x), 
y = f(w) = f(g(x)). Then we can write 


 — © (F(g(a))) = F(o(2))o'(@) = SE. 


Machine: 


we have gained 3 basic rules of differentiation, 


1. sum rule: 


d df — dg 
+ (F(a) +(e) = £4 F 
2. product rule: E F r 
© (F(aala)) = 9a) E+ Fay, 
3. chain rule: q df d 
= g 
5 lala) = FF. 
Extras: . 
a) = px. 


2.2 Differentiation of trigonometric functions 


How do we find the derivative of f(x) = sina? By definition 2.1 (on pg. 22), the derivative 
of f(x) = sin at the point x = c is 


fy = tim fe+N= 10 
dx h-0 h 
ay sin(c + h) — sin(c) 
h-0 h 
Bs Spe sin(c) cos(h) + sin(h) cos(c) — sin(c) 
h-0 h 
‘ : cos(h) — 1 sin(h) 
= lim {sin(o h + cos(c) 1 \ 
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If we know the limits of sin(h)/h and (cos(h) — 1)/h as h > 0, then we will know f’(c), 
since c is unrelated to h. 


First let us consider sin(h)/h. We draw a circular sector with a very small angle, where 
the curved side is of length «x (in radians). For small angle x (in radians), sinz and x are 


almost equal, i.e. 


sin x 
xl. 
x 
1 
x 
|| 
1 
Figure 2.12: Very small circular sector. Right angled triangle with hypotenuse and 
adjacent approximately equal t.e. two radii from circular sector. 


As h > 0, sin(h)/h — 1, if you work in radians. Actually 


ee sin(h) fe sin(0 + h) — sin(0) 
ho A h0 h 
= —(sinz) 
d «=0 
= cos(0) 


that is, the derivative of sinxz at x = 0 is one, or the tangent line of y = sinx at x = 0 is 


y=. 


(a) Tangent sinx at x =0. (b) Tangent cosx at x = 0. 


Figure 2.13: Examining the tangents at x =0 gives an insight of the derivative at this 
point. 
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Second, we notice that 


kim cos(h) — 1 od. TR, cos(0 + h) — cos(0) 
h—0 h h—0 h 
d 
= ag (68 L) Le 
= sin(0) 


i.e. the derivative of cosx at x = 0 is zero, since the tangent line of y = cosx at x = 0 is 
horizontal (y = 1). So 


kim cos(h) — 1 =, 
h->0 h : 
Therefore 
f"(e) = cos(c), 
ie., 


Similarly, we can derive 


ag (8 £2) = —sing. 
NOTE: The change in sign when differentiating cos z! 


Exercise 2.1. Now that we know the derivatives of sin and cos, we should be able to 
calculate the derivative of tana. Try it yourself before next lecture!* 


“End Lecture 8. 
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Since we know the derivatives of sin and cos, we should be able to find the derivative of 
tan since : 
sin x 
tanz = 


COs x" 


Before we attempt to differentiate (using the product rule), there is another “so called” 
rule called the quotient rule, really, it is simply the product rule for rational functions. We 
derive the formula below, however if you do not remember the formula, it’s best to work 
with the product rule. 


First we write 


as a product, then apply the product rule as follows: 


“ ae) = f'(x) [g(x)J* + f(a) {[a(a)1""} 


We need the derivative of [g(x)|~'. This is obtained by the chain rule, if we consider 
r(u) =1/u, then r’(u) = —1/u?, so 


| 
ey 
Pe 
= 
2 

I 

| 
Zz 
Se 
= 
S 
hoo 


Finally, we have 


d 9) —praiencnn ey GO) 
= [43] = r@@r - Fe) 


finding a common denominator, we may write 


d [f(x)) _ f'w)g(@) — f(@)9'(@) 
mabey oor er 
This is known as the quotient rule. 


Example 2.18. Consider the function tanz. We will use the quotient rule where f(x) = 
sinx and g(x) = cos, therefore f’(x) = cosx and g'(x) = — sina. Hence 


re ao (=) 
dx dx \cosx 
cos? x + sin? x 
cos? x 


1 
cos? x 


= sec? x 


= 1+tan?z. 


Exercise 2.2. Try find the derivatives of cosecx, secx and cot x. 
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Example 2.19. Consider the function cos(x? + 2x7). Let us choose f(g) = cos(g) and 
g(x) = x? — 22, i.e. f'(g) = —sin(g) and g'(x) = 2x — 2. Then, using the chain rule 


d d 
re [cos(ax? + 22) == Ns [f(9(x))] 
= f'(g(a))g'(2) 


=(Q2= 2). (sin(a? = 22), 


4 


Example 2.20. Consider the function cos* x sinz. Here we will first apply the product 


rule, so that 
d d 
7 [cos* x sina] = sine (cos! x) + cos4 (sin az) 
x 4s x 
The derivative of sinx is easy, we know this. All we need to do now is apply the chain 
rule on cos*x. Here we will choose f(u) = u*, therefore f’(u) = 4u3 and g(x) = cosa, 


whose derivative is g/(x) = —sinx. Thus we have 
Ee [cos* x sin x] = EPR Ceny) + cos’ x - cos x 
dx dx 


= sina(f"(g(x))9!(w) + 608° « 
= sinx-4cos’ x: (—sinx) + cos’ x 


= —Asin? xcos* x + cos” x. 


2.2.1 The derivatives of inverse trig functions 


Now that we know the derivatives of the trigonometric functions, let us consider the 
inverses of these functions, i.e. sin7! a, cos! a and tan~! x. 


NOTATION: sin~! & is the inverse function of sina, whereas (sinz)~! = 1/sinz. 


First let us see how to define sin~! x. To define sin~!, we first must put a constraint on the 


domain of sin. Recall, a function which associates more than value of x for f(x) cannot 
have an inverse. So we choose the domain [—4, 5]. 


1 


sin“! : [-1,1] > [—, 3]. 
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Let f(u) =sinu and g(x) = sin“! z, we need to work out g'(x). We know that f(g(x)) = 2, 
so differentiating on both sides we have 


, ! = (a) = : T= : 
fe@w'@=1 = I@)=aEy 


cos(sin~! x)” 


If we let 6 = sin7! x, then we have sin@ = x. Now we construct a right angle triangle 
which satisfies sin? = x. 


V1— 2? 


x 


Figure 2.15: If sin(@) = x, one can construct such a right angle triangle using 
Pythagorus’ theorem. 


Then we know from the triangle that 


cos@ = V1— 22. 


Therefore 
5, (sin xr) = 


In the same way, we can work out 


Exercise 2.3. Define g(x) = cos”! x first and then work out g'(z).° 


5Fnd Lecture 9. 
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2.3. Maxima, minima and second derivatives 


Consider the following question: given some function f, where does it achieve its maximum 
or minimum values? 


First let us examine in more detail what f’(x) tells us about f(z). 


a 
f'>0 f' <0 f'>0 f' <0 


a) Maxima and minima of some graph, y = f(x). 


(b) Positive slope. (c) Negative slope. 


Figure 2.16: Displaying the change in first derivative before, at and after maxima and 
minima. 


We notice that: 


If f(x) is increasing, then f’(x) > 0, 


If f(a) is decreasing, then f’(x) < 0. 


Therefore, troughs and humps occur at places through which f’ changes sign, i.e. when 
f’ = 0, where 


Trough = local minimum, Hump = local maximum. 
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The derivative gives us a way of finding troughs and humps, and so provides good places 
to look for maximum and minimum values of a function. 


Example 2.21. Find the maximum and minimum values of the function 


3 3 
f(x) = x — 3x, onthe domain — 5 <a< 5 


Differentiating f(x) we have 


f(x) = 327 — 3 = 3(a7 — 1) = 3(a + 1) (x — 1). 


Hence, at x = +1, we have f’(x) = 0. The value of the function at these points are 


It also maybe possible for the function f to achieve its maximum or minimum at the ends 
of the domain. Therefore we calculate 


So comparing the four values of f, we know where f(x) is biggest or smallest. So we 
finally conclude 
f(1) =—-2 (min), f(—1)=2 (max), 
i.e. we have a minimum at the point (1,—2) and a maximum at (—1,2). 
If we had been looking in the range —3 < x < 3, then at the ends 
f(3) =18 (max), f(—3) =—18 (min). 


To draw the graph of f, we need to find the values of f at some important points, such 
as when « = 0, when f(x) = 0 and f’(a) = 0. We already know the points at which 


f(x) = 0. 
When x = 0 we have f(0) = 0, ie. the graph passes through the point (0,0). 
When f(x) = 0, we need to solve the equation 

2° — 32 = 2(x* — 32) = 0. 


So either 7 = 0 (we already have this point), or 


x —3=0 = t= +V3. 
So the graph also passes through the points (0, V3) and (0, —V3). 
Now let us examine the derivative further, f’(x) = 32? — 3. 


f(z)>0 when x<-l, 
f'(z) <0 — when -l<a2<l, 
f'(cz)>0 when e>1. 
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Figure 2.17: Sketch of f(x) = x? — 3x, displaying turning points (maximum and 
minimum). 


Sketching graphs, things to remember: 


1. Find f(x) when x = 0, i.e. where the graph cuts the y-axis. 
2. Find « when f(x) = 0, i.e. when the graphs cuts the z-axis. 


3. Find x when f’(x) = 0, ie. the stationary points of the graph. (Plug into f(x) 
to find y value). 


4. Determine the sign of f’(a) on either side of the stationary points to determine 
weather stationary points are minimum, maximum or points of inflection. 


Basic principles: 


Let f : [a,b] > R. If f achieves a local maximum or local minimum at x, then either: 


(i) f(x) =0, or 
(ii) «=a, or 
(iii) « =b, or 

) 


(iv 


So to find the local maximum/minimum of f, it suffices to list possibilities in (i)-(iv) and 
then check. 


Example 2.22. Consider f(x) = |x|, on the domain —1 < a < 1. We know f’(0) doesn’t 
exist, but at « =0, f(x) achieves its minimum value of 0. 


Example 2.23. Consider f(x) = 1/x on the domain —2 < x < 2. There is no maximum 
or minimum on this range. f’(2) = —1/2? is not defined at 2 = 0 (it actually tends to 
too either side of x = 0). 
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Example 2.24. Consider f(x) = x° on the domain —2 < x < 2. f’(x) = 3x”, which is 
zero if x = 0. So zero is a stationary point for x°, but it is neither a hump or a trough. 
It’s a point of inflection. 


y=-# y=u 


By 


ee, A 


(a) y = || (Ex. 2.22). (b) y = + (Ex. 2.23). (c) y= a? (Ex. 2.24). 


Figure 2.18: Different options for when f'(x) = 0. 


2.3.1 Second derivative 


To characterise troughs and humps (local maximums and minimums), we need the knowl- 
edge of second derivatives. 


If we start with some function, say f(x) = x?, we know the derivative f’(a) is well defined 
for every x, and f’(x) = 2a. We can view f’() itself as a function, so we may differentiate 
it again to have 
de> 9) d? 7 
£(¢(@) = SU@) = F'@) =2 


In this way, we can define f”, f®), f“ and so on. 


Example 2.25. 


f(x) = x", n is a positive whole number. 

fie) = na", f(a) = n(n — 1a", 

f(x) = n(n — 2)(n — 2)2"-3, ... 

f(x) = n(n —1)(n—2)...2-1=nl, recall ! = Factorial. 


f(x) =0, m>n. 


In physics, if f(¢) represents distance as a function of time, then f’(t) represents speed 
(i.e. the rate of change in distance), and f”(t) represents acceleration, i.e. the rate 
of change of speed. This is the key to understanding how things move, in particular 
under gravity. 
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The geometric interpretation of f”: 


1. If f” > 0, then the slope of the tangent line is increasing in value (from left to right), 
so possible shapes for f(a) are like: 


Figure 2.19: Slope of tangent increases in value (possibly negative to positive). 


Therefore if f’(c) = 0 and f”(c) > 0, then around c, f(x) is a trough and we can 
expect a local minimum value of f at x =c. 


Example 2.26. 
f(z) = 2°—3z, 
f@) = 327 =3858@=1)=32@-1)e+h, 
Fey Si.6z: 


At x =1, f’(x) =0, f”(x) = 6 > 0, so we have a trough at x = 1. 


2. If f” < 0, then the slope of the tangent line is decreasing (from left to right), so 
possible shapes for f(a) are like: 


Figure 2.20: Slope of tangent decreases in value (possibly positive to negative) . 


Therefore if f’(c) =0 and f”(c) < 0, then around c, f(a) is a hump, and we expect 
a local maximum value of f at c. 


Example 2.27. Continuing with f(x) = #3 — 3z, at = —1 we have f’(—1) = 0, 
f"(-1) = -6 <0. 
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Figure 2.21: Characterising the turning points using second derivatives of 
— rd 
y= x? — 3x. 


3. If f’(c) = 0 and f”(c) = 0, then the slope often doesn’t change sign, i.e. it goes 
from positive slope to zero to positive slope (decreasing to zero then increasing), or 
negative to zero to negative (increasing to zero then decreasing). These points are 
points of inflection and possible shapes are like: 


c 


Figure 2.22: Points of inflection at c, slope goes from positive to positive or negative 
to negative. 


Example 2.28. f(x) = x, so f’(z) = 327 and f(x) = 6x. At « = 0 we have 
#0) = f"(0) = 0. 


WARNING: However, f(x) doesn’t necessarily need to be a point of inflection when 
the second derivative is zero. 


Example 2.29. Consider f(x) = x*. Now we have f’(x) = 4x3 and f"(x) = 122°. 
Notice that f’(0) = f’(0) = 0, however when we sketch f(a), we realise it has a 
local minimum at x = 0. 
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—_> 


x 


Figure 2.23: Graph of y = x*, clearly doesn’t have an inflection point at x = 0, yet 
f'"(0) = 0. 


Therefore we say that if f has a point of inflection at « = c, then f’(c) = f”(c) =0 
is a necessary condition. 


Example 2.30. (A slightly more advanced example). Consider the function 


a—4  (4—2)(x+2) 


NO et Galant) 


So f(x) = 0 when the numerator of f(x) is equal to zero. That is 


(a — 2)(x+ 2) =0 v=2,"=-2. 


Something strange happens at the points « = 1 and x = —1. We get zero in the denomi- 
nator and the numerator is negative in both cases so f(+1) + —oo. 


The graph cuts the y-axis when x = 0, that is 


Oe 


=—=4. 
—*l 


Now let us calculate the derivative, it is a rational function so we may use the quotient 
rule (or apply the product rule), so we have 


1, _ 2a(@?-1)-2x(a?-4) 6 
f(x) (x2 a 12 (x? = 1)?’ 


So f’(z) = 0 when the numerator is zero i.e. when 62 = 0, thus there is only one turning 
point at x = 0. 


Now let us differentiate f’(x), again using the quotient rule we have 


_ 6(a? — 1)? — 6x(2- (a2 -1)- 2x) — 6(a? — 1)? — 242? - (x? - 1) 


i 
ro @ i) @ 1) 
It remains to check the value of f” at the turning point « = 0. Therefore 
6(—1)? —0 
"(0) = ——,__ = 6 > 0 
f"(0) (<ijt > 0, 


and we have a local minimum at this point. 
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Figure 2.24: Graph of y = (a? — 4)/(a? — 1), i.e. a rational function with only one 
turning point (minimum) at x = 0. 


At x = +1, we have what is known as asymptotes, lines the graph never touches, but gets 
very close to as we approach plus or minus infinity.® 


®Fnd Lecture 10. 


CHAPTER 2. DIFFERENTIATION AT 


2.4 Real life examples 


Example 2.31. You are required to design an open box with a square base and a total 
volume of 4m? using the least amount of materials. 


We are hunting for particular dimensions of the box. Lets give them names, so we will 
call the length of the base edges a and let the hight of the box be called b. See Fig. 2.25. 


a 


Figure 2.25: Square based open top box. 


In terms of a and b, the volume of the box is 


a’b = 4 (m’). 


We want to minimise the amount of material, i.e. the question is how much do we use? 
Consider what we know about the box. We have a base of area a? and 4 sides, each with 
area ab. So the total area is 

a? + dab. 


We want to find the minimum value of a? + 4ab, subject to the condition a?7b = 4 and 
more obviously a > 0, b > 0. From the condition, we see that b = 4/a?, so we can rewrite 
the amount of material in terms of a: 


4 1 
+ dab =o + 4a) =o + S 
a 


So now we can write a function for the material area solely in terms of one of the lengths 
of the box (the other is now fixed by using the condition). That is 


Now written like a function, it is easy to see how we would minimise the material, that is 
by minimising the function f(a). So we have to calculate f’(a), which is 


16 
f'(a) = 2a — G2" 
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and now we simply need to see when f’(a) = 0. Finally, we have 


16 


2a-—=0 = > 2° =16 a= 2 
a 


The only real number at which f’(a) = 0 is a= 2. The function f’(a) makes sense except 
at a = 0, which is outside the range (since a > 0). In addition, in this situation we don’t 
have a closed domain (i.e. no ends) since 0 < a < +00. We can differentiate f’(a), which 
gives ay 

f"(a) = 2+ aa 
This is positive for the whole range (and at a = 2), along with sketching the graph in Fig. 
2.26, we can mathematically confirm choosing such an a would use minimal materials. 


Figure 2.26: Graph of f(a) =a? +16/a, displaying changes in area of material 
required for different choices of length a, given fixed volume condition on b. 


NOTATION: For an open domain, for example a < x < b, we use curly brackets, i.e. 
a € (a,b). For an closed domain, for example a < x < b, we use square brackets, i.e. 
a € [a, b). 


Example 2.32. We have a pair of islands, island 1 and island 2, 20Km and 10Km away 
from a straight shore, respectively. The perpendiculars from the islands to the shore are 
30Km apart (along the shore). What is the quickest way between the two islands that 
goes via the shore? 


We are trying to find a point along the shore, which we want to visit when going from 
one island to the other. This point can specified by the distance from the perpendicular 
of island 1, call this distance x. 
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yA 
eisland 1 
PS 10? + (30 — x)? 
ey island 2 
20 ae Pai 
is 202 + 22 
r baa y a x 
x 30-2 


Figure 2.27: The set up of islands relative to the shore (x-azis). 


The total distance D(x) to be travelled is given by the formula 


D(x) = V 20? + «? + +/10? + (30 — z)?. 


On geometric grounds, we see that 0 < x < 30, and D(z) is differentiable in this range. 


1 
D(z) = 5(20? + 22) ban 5 (10? + (30 — x)?] 2 2(30 — z)(—1) 
x 30-2 
V20% +22 /(30—2)? +102 


To find the minimum distance we set D’(x) = 0, so we have 


x - 30-2 
V202 + a? (30 — x)? +102 


Squaring both sides we get 


a (30-2)? 
202 +a2 (0-2)? +102 


(x”)[(30 — x)? + 107] = (30 — 2)?(20? + x?) 
Expanding the brackets we see 

102? = 207(30— 2)’. 
Take the square root of both sides of the last equation gives 


10z2 = +20(30 = x). 


Since 30 — x > 0, “-” sign is not possible, so 
102 = 20(30-—2) => «=20. 


The only possible places where minimum can occur are 


where 


D(0) = 20 + /10? + 30? = 51.6, 
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D(20) = 20V2 + 10V2 = 30V2 = 42.4, 


D(30) = 20? + 302 + 10 = 46.06. 


So the minimum value does occur at x = 20. 


yh 


qisland 1 
island 2 
20 ? 


1 aes we ‘. 


20 10 


Figure 2.28: Optimum route from island 1 to island 2, arriving and leaving shore at 
an angle of 1/4 radians. 


The optimal route is to leave the shore at the same angle of arrival. 


yh 


gisland 1 


20 5 ai 2 


Ry 


Figure 2.29: Sketch of alternative method of calculating shortest distance between the 
islands, using simple geometric properties, gaining same result. 


Could we have deduced that this route was the shortest without calculus? Yes! We could 
have reflected island 2 in the shore line to obtain an imaginary island. Then it is easy to 
see that the shortest route from island 1 to the imaginary island is a straight line.” 


"End Lecture 11. 


Chapter 3 


Exponentials and Logarithms 


3.1 Exponentials 


An exponential function is a function of the form 


f(x) =a", 

where a is a positive constant. 
Example 3.1. 

yt 

1 a=2>1 

ft | 
z 
(a) y = 2”. 
Figure 3.1: The domain: —oo < x < +00; the range: (0, +00). 


For 


a> 1, f(x) increases as x increases. 


a <1, f(x) decreases as x increases. 


a=1, f(x) =1. 


a° = 1 for each a, so the graph always passes through the point (0, 1). 


ol 
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3.1.1 Slope of exponentials 


First let us consider the slope at x = 0. 


Example 3.2. Suppose we have f(x) = 2”, then applying the definition of the derivative 
we have 


‘ fO+h)—fO)_ |. 2- 
pS h ah ah 
h | f(s 
0.1 0.7177 
0.01 0.6955 
0.001 0.6933 
0.0001 | 0.6932 


So for f(x) = 2”, (a = 2), we have slope ~ 0.693 at x = 0. 


Similarly, for f(x) = 3”, (a = 3), we have slope ~ 1.698 at x = 0. 


Therefore, we expect that there is a number between 2 and 3 such that the slope at 
x =0 is 1. This number is called e, where e ® 2.718281828459.... The number e is 
irrational. 


Figure 3.2: Graph of y =e”, which has tangent with slope of 1 at x =1. 


The fact that the slope is 1 at « = 0 tells us that 


Therefore, if h is small, then e” — 1% h, ice. 


e* 1+h. 
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We call f(x) = e® = exp(x) the exponential function. 


To find the slope at x = c, we need to look at 


ye, fetA— FO) eth eet), eh, 
ce ne a ar = a ae 
i.e. the derivative of e* is itself, 
d xv x 
an )=e . 


Example 3.3. Consider f(x) = eY!+*. Here we will employ the chain rule. Choose 
g(x) = /1 +2 and f(u) =e", so we have g/(x) = $(1+ x)-3 and f’(u) =e". 


f(T) = f(g(2))a!o) 


1 
e meat + a2 


e 1+a 


2W1l+ta 
3.2. The natural logarithm 


Consider the inverse function of f(x) = e*. In f(x), every positive number occurs as the 
exponential of something, i.e. M = e! for an appropriate t. The number ¢ is called In(M): 
the natural logarithm of M. In other words 


“In M is the number whose exponential is M”: e™™ = M. 


In this way, we define the function of the natural logarithm 


g(M) = 1nM. 
Example 3.4. 
yh 
t z 
(a)y=e* (b) y=Inz 
Figure 3.3: For nx the domain: (0,+00); the range: (—o0, +00). The graph y = e” 
has a horizontal asymptote at y= 0, while y = In has a vertical asymptote at x = 0. 
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By definition, 


e™M™ _ M, In(e*) =a, 


which means that if you perform In(exp) or take the exp(In), then we get back to where 
we started. 


3.2.1 Characteristic properties 


1. n(MN)=InM4+InN. 
2. In(M?) = piInM. 


Logarithms are used among other things to solve “exponential equations”. 
Example 3.5. Find x, given 3” = 7. Taking the logarithm of both sides we have 
In(3*)=In7 = gln3=1n7. 


Rearranging we have 


oa 8 1 i 
In 1.1 
Exercise 3.1. Show that 
ais r)= 
dx 7 


Hint: put y = Inz. 
Example 3.6. Consider the function f(x) = In(cosa). Choose g(x) = cosa and f(u) = 


Inu, so we have g/(a) = —sinz and f’(u) =1/u. Thus 
d 
qalin(cosx)) = fi(9(a))9"(@) 
by 
i 
Se (— sin x) 
= —tanz. 


Similarly we have 


dd, 1 
q, Sin(n 2) = cos(In x) - = 


Differentiation of other exponentials 


In order to differentiate for example 3”, we must express it in firms of e” : 
3= eins — 3% (Cama = erins 


Therefore we calculate the derivative of 3” as follows: 


d x = d aln 
pe gee? 
f'(g(2))g'(2) 
= erln3 . In3 


= 3*-In3. 
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Here we chose g(x) = x1n3 and f(u) = e” so that g'(x) = 1n3 and f’(u) = e”. 
In general, for any positive constant a 


d 


aa? =a" Ina. 


NOTATION: log, x = Ing, is the “proper” way of writing the natural logarithm. 


3.2.2 Logarithms base a 


We can also define log, (x) to be the number m, i.e. log,(%) = m is such that a” = x. In 
this way we can think of logarithms as a different form of writing powers. 


Example 3.7. log;)(1000) = 3. 


The derivative of log, (x) is 


© (Woga(2)) = 


x In(a) 


Exercise 3.2. Try to show the above statement is true. Hint: use the chain rule.! 


1nd Lecture 12. 
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3.3. Exponential growth and decay 


Let y = f(t) represent some physical quantity, such as the volume of a substance, the 
population of a certain species or the mass of a decaying radioactive substance. We want 
to measure the growth or decay of f(t). 


Definition 3.1. We define the “relative rate of growth (or decay)” as 


actual rate of growth (or decay) _ f'(t) _ dy/dt 


ete 0 ae CH 


In many applications, the growth (or decay) rate of the given physical quantity is constant, 


that is 
dy/dt dy 
— =a, or —=ay, a=constant. 
y dt 


This is a differential equation whose solution is 
y(t) = ce™, 


where constant c is determined by an initial condition, say, y(0) = yo (given). Therefore 
we have 


y(t) = yoe™. 
This means that if you start with yo, after time t you have y(t). 


If a > 0, the quantity is increasing (growth), 


If a < 0, the quantity is decreasing (decay). 


3.3.1 Radioactive decay 


Atoms of elements which have the same number of protons but differing numbers of neu- 
trons are referred to as isotopes of each other. Radioisotopes are isotopes that decompose 
and in doing so emit harmful particles and/or radiation. 


It has been found experimentally that the atomic nuclei of so-called radioactive elements 
spontaneously decay. They do it at a characteristic rate. 


If we start with an amount Mo of an element with decay rate 4 (where A > 0), then after 
time t, the amount remaining is 
M = Moe. 


This is the radioactive decay equation. The proportion left after time t is 


M 
Be gates 
Mo 
and the proportion decayed is 
M 
1-—- —=1- ea, 
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3.3.2 Carbon dating 


Carbon dating is a technique used by archeologists and others who want to estimate 
the age of certain artefacts and fossils they uncover. The technique is based on certain 
properties of the carbon atom. 


In its natural state, the nucleus of the carbon atom C!? has 6 protons and 6 neutrons. 
The isotope carbon-14, C', has 6 protons and 8 neutrons and is radioactive. It decays 
by beta emission. 


Living plants and animals do not distinguish between Cl? and C™, so at the time of death, 
the ratio C!? to C™ in an organism is the same as the ratio in the atmosphere. However, 
this ratio changes after death, since C™ is converted into C!? but no further C' is taken 
in. 


Example 3.8. Half-lives: how long before half of what you start with has decayed? When 
do we get M = 5Mo? We need to solve 


M 2° ’ 


taking logarithms of both sides gives 


n(3) = M t = 


So, the half-life, T)/2 is given by 


In (2 
Lg = ‘ y 


If A is in “per year”, then 7/2 is in years. 


Example 3.9. Carbon-14 (C1!) exists in plants and animals, and is used to estimate 
the age of certain fossils uncovered. It is also used to trade metabolic pathways. C4 
is radioactive and has a decay rate of X = 0.000125 (per year). So we can calculate its 
half-life as 


In2 
Ty jo = ——— © 5545 years. 
/2 ~ 0.000125 ee 


Example 3.10. A certain element has T}/2 of 10° years 


1. What is the decay rate? 


_ In2 _ 0.693 


N= aw 
Tij2 108 


= 7x 10-7 (per year). 


2. How much of this will have decayed after 1000 years? The proportion remaining is 
M a pAb ag TROT R10? _. 7X10 9993. 
Mo 


The proportion decayed is 


M 
1- — 1-0. = 0. : 
Mo 0.9993 = 0.0007 
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3. How long before 95% has decayed? 


M -7 
= 1—0.95 = 0.05 = —7x10~'t 
Mo 0.95 = 0.05 = e ; 


taking logarithms of both sides we have 
In(0.05) = —7 x 107"t 
which implies 


In(0.05) _ —2.996 


fa 6 
SEeIOne eee 4.3 x 10° (years). 


WARNING: Half-life Ti/2 of a particular element does not mean that in 2 x Tj/2, the 
element will completely decay.” 


?End Lecture 13. 
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3.3.3 Population growth 


Example 3.11. Suppose a certain bacterium divides each hour. Each hour the population 
doubles: 


4 5 
16 | 32 


| Hours 1] 2 
| Population 2/4 


After t hours you have 2° times more bacteria than what you started with. In general, we 
write it as an exponential form. 


If a population, initially Py grows exponentially with growth rate \ (where A > 0), then 
at time t, the population is 
P(t) = Poe. 


Example 3.12. Bacterium divides every hour. 


1. What is the growth rate? 


We know that when t = 1 hour, we are supposed to have 
P=2Pp, 


sO 


2Po = Poe™? oe \ = In2 © 0.693. 
2. How long for 1 bacterium to become 1 billion? 
Py His X= 0.698, P10", 
therefore we may write 
P=Poe™ => 109 = ¢0-698¢, 


taking logarithms of both sides and re-arranging for t, we have 


= 91n 10 


~ h ; 
0.693 30 hours 


3.3.4 Interest rate 


An annual interest rate of 5% tells you that £100 investment at the start of the year grows 
to £105. Each subsequent year you leave your investment, it will be multiplied by the 
factor 1.05. 


In general, if you initially invest Mp (amount) with an annual interest rate r (given as 
percentage/100), then after t years you have 


A= Mo(1 + r)', 
where A is the future value. We could write this as an exponential as follows: 


A(t) = Moe** = Mo(1+ fi 
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Taking logarithms we have 
At =tin(1+r), 


so we may write 
A(t) a Maelo, 


ASIDE: In fact, for small r, In(1 +r) r. 3 


°Fnd Lecture 14. 
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Chapter 4 


Integration 


4.1 The basic idea 


We are interested in calculating areas under curves. 


Example 4.1. 


a b £ 


(a) Shaded area under y = f(x). (b) Rectangular strips under y = f(x). 


Figure 4.1: Integration calculates the shaded area under the curve y = f(x), we could 
do this by dividing the area into rectangles and summing the area of all rectangles. 


How do we do it? 


1. We divide the interval a < x < b into pieces (say equal length). 
2. We build a rectangle on each piece, where the top touches the curve. 
3. We calculate the total area of the rectangles. 


We say, if the division is very fine, we will get a good measure of the area we want. We 
watch what happens as we make the division of the “strips” finer and finer. 
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Example 4.2. Consider the function f(x) = x on the interval 0 < #7 <1. 


i 
n 


Figure 4.2: Integrating under the curve y= 2x, from x =0 tox =1. 


We divide [0, 1] into n equal pieces. The divisions occur at 


1 2 k-1k n—-1 m4 
nn?” n eh n on * 


0 


We have n+ 1 points and we put a rectangle on each point. The rectangle between nt 


and © will have height f(*) = © (see Fig. 4.2), and the area of this rectangle is 
k 1, 
te on 


BY ~Y 
height width 


The sum of the area of all rectangles on the interval is 


fi 2 k n 1 

ge eb eg Si ia: en tk 4 + n) 
-- eae 
— ne 2 


I 
N] rR 
a 
me 
+ 
Sie 
SS 


and 
1 1 1 


~{1+-—- ~ : 
s(1+2) 45, as n > CO 


1 


That is the sum is approaching the actual area 3. 


Therefore, as we increase n so as to get finer divisions, the area approaches the exact area 
under the curve. You could do this for all functions, however, there is a much quicker way. 


Idea: we want to find the area under the curve, let us call it y = f(t). 
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Example 4.3. 


a b t a x b t 
(a) Shaded area under y = f(t). (b) Integrating y = f(t) up to point x. 


Figure 4.3: Writing the shaded area under the curve y = f(t) as a function of x, where 
x is point within the desired interval [a, b]. 


We think of the area as a function of x, say A(x). If we know A(x), then we know the 
area, i.e. A(b) — A(a). So we want to find this function A(x). 


We don’t know A, but we can say something about it. Think about how it is related to 
f. What we know is that 


A'(x) = f(a). (4.1) 
Why can we say this? We need to understand what happens to 
A(a +h) — A(x) 
h 


as h > 0. (4.2) 


t=c t=a2+h t 


Figure 4.4: Consider one “strip” under the curve y = f(t), to find the limit (4.2). 


The difference A(x + h) — A(x) is the area between t = x and t = x +h. So the area is 
roughly rectangular (if h is small) with height f(a) and base h. So the area is approxi- 
mately f(x) -h. Therefore 


A(a+h)—A(a) = f(z)-h => Aenea) ee f(x), 
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and 


A(xz+h) — A(x) 
h 


64 


> f(x) ash 0. 


Thus, by the definition of the derivative, we have A’(x) = 
the antiderivative of f(x).+ 


f(x). We defined A(z) as 


1nd Lecture 15. 
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Notice that if A(a) is an antiderivative of f(x), then A(x) + C, where C is any constant, 
is also an antiderivative of f(a). 


We employ the following notation: 
[f@de= Fe) +e, P(e) =s@). 


J f(x) dz is referred to as an indefinite integral, where f(x) is called the integrand. 


4.2 Fundemental Theorem of Calculus 


Example 4.4. Consider f(t) = t?. We want to calculate the area under y = t? between 
O and 1. 


> 
t 


0 x I 


Figure 4.5: Integrating to find the shaded area under the curve y = t?. 


The area up to x is represented by 


A(a) = t dt. 


We know A’(x) = x? = f(x). We need to find A(x). Can we think of a function whose 


derivative is x7? Yes! Consider ae In fact, for any constant C, 


d (1. 9 


Our function A(x) is one of the family of functions 50° +C, but which one? 


We need to fix the constant C. We have A(0) = 0, so if A(z) = $23 + C then we must 
have C' = 0. So finially A(x) = 52°. Hence, A(1) = $, which is the answer to the original 
question. 


We can interpret the above result as 


o # dt = A(1) — A(0), -A'(z) = 2. 
0 
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This is actually the fundamental theorem of calculus, which says: If F’(x) = f(x) between 
a and b, then 


b 
/ f(x) dx = F(b) — F(a). (4.3) 


The integral ft f(x) dz is called the definite integral of f(x) from a (the lower limit) to b 
(the upper limit). 


ie f(x) dz is the limiting value of the sequences derived from the method of divisions 
when calculating the area under the curve, described previously. Thus, the integral may 
represent areas under curves. 


Example 4.5. Suppose we want to integrate the function x* over the interval [2,3]. That 
is, we want to calculate 
3 
| a! dx. 
2 


Think of a function whose derivative is z+. The answer is a So, we write 


3 3 
ah adr = Ls = bs 195 — ay 
2 5 9 5 5 5 


Example 4.6. Suppose we want to integrate the function 1/x? over the interval [1,2]. 
That is, we want to calculate 


If we put F(x) = —1/z, then 


So we can write 


ie 1 1]? 1 yt dl 
5 dx = = =-. 
1 2 Ly 2 1 2 
The integral { i + dx represents the area under the curve y = + between 1 and 2, therefore 
we understand that this integral makes some geometrical sense. 


Figure 4.6: Integrating to find the shaded area under the curve y = + on the interval 
[1,2]. 
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Example 4.7. Consider the integral of 1/x?, but this time on the interval [—1, 1]. 


Figure 4.7: Integrating to find the shaded area under the curve y = + on the interval 
[—1,1]. However, the curve has a vertical asymptote at x = 0. 


However, the area under the curve in this interval is not —2! What is wrong here? Think 
about this for next time.? 


?End Lecture 16. 
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Previously, we chose an antiderivative which is correct for the given integrand 1/x?. How- 


ever, recall 
d 1 to 


That is F’(a) = f(a) doesn’t hold for —1 < a < 1. We have to be sure the function is 
well defined over the entire interval over which we integrate. 


Example 4.8. Consider the function f(x) = 2°. 


[—1, 1] is 
1 
[ ea- Fa an 1 6, 
1 4 1, 4 4 


Then the integral over the interval 


po Pee 


Figure 4.8: Integrating to find the shaded area under the curve y = x® on the interval 
[—1,1]. Have two shaded regions bounded by the curve and the x-axis. 


In this case, the area cancels out. The shaded area is actually given by 
1 0 7 1 1 1 0 
if a? dx + / x? dx| = a si fa 
0 -1 a ge ee shat 


4.3 Indefinite integrals 


ere 
"4° 4 


So far we know 
[to dp= Fig) 4C, F'(z) =F (@): 


What function F can we differentiate to get f? 


Powers of 2: 


n = 1 n+1 oe 
fe aa +C, (n#-1), 
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d grtl a 
= (=) aide 


1 
[gaeamlelte, 
x 


since 


For n= —-1, 


since for x > 0 


d ik 
qq ln) = ey 


and for x < 0, 


Special rule: 
Let us consider the derivative of the logarithm of some general function f(x), i.e. 


d 1 d _ f(z) 


An [In(f(x))] = fx) ae [f(x)] f(a) 


This implies that 


i sae 
[aga mue@sre 


where C is some arbitrary constant of integration. Hence, if we can calculate integrals by 
inspection if the integrand takes the form f’(x)/f(ax) by the above formula. 


Example 4.9. Consider the the following integral: 


27 +5 
PS ae 
i e+ 5e +3 o 
Now, if we choose f(x) = 2? +5a +3, then f’(x) = 27+5. So, if we differentiate In(f(z)), 


in this case we have 
24 +5 


x2 +52 +3’ 
by the chain rule. Thus, we know the integral must be 


d 
— [In (a? + 5a +3)] = 
xv 
ES In(e? + 52+ 3) EC, 
where C’ is some arbitrary constant of integration. 


Example 4.10. Consider the following integral: 


I= | ——dz. 
are hee 

Now, if we choose f(a) = 2a+2 then f’(a) = 2. However, the numerator of the integrand 

is 3. Not to worry, as we can simply re-write or manipulate the initial integral as follows: 


3 1 2 3 2 
i= 8 j= ae 
las ‘i fas is sla e 


Since 3/2 is a constant, which we are able to take out of the integral sign, we need not 
worry about this and can proceed with the integration using what we have learnt above, 
giving 
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To check, we differentiate the above expression, so 


dI d |3 3 1 
ae ae 5 Ina t2)+C 


which is correct! 


This “special case” is an example of a method called substitution, and is not limited 
to integrals which give you logarithms. Nevertheless, it is a good sighter for what’s to 
follow. Together with inspection, it can be extended for other function by choosing suitable 
substitutions (i.e. f(x)). 


Trigonometric functions: 


d 
[cos dx =sinx+C,_ since 7 sin x) = cosa, 
L 


d 
[sucae =—cosx+C, since aa cos x) = sina. 
2 


Exponential function: 


ie dzx=e*+C, since —(e*) = e”*. 
x 


As with differentiation, we also have a sum rule for integration, that is 


[ro + g(x)) dx = [ito dx + [oe dx. (4.4) 


In other words, the antiderivative of a sum is the sum of the antiderivatives. 


Also, 
[ Kt) dx = K | He) dx, (4.5) 


where K is a constant, i.e. a constant can be taken outside of the integration sign. 


Example 4.11. 


[ext +2045) de = [ostact [rede + [sac 
3 fatar+2 faders f dz 


NOTE: Do not forget the constant C when the integral is indefinite! 


Example 4.12. 


[@-ve@t+ana = [oe-ot+ 20 - 2a 


[ood fotar+2 fade— f rae 
1 
7 


1 2. 
gt a | 30 2e +C. 
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Example 4.13. 


vt +1 a 1 
i 2 dz = [G | =) dx 
il 


4.3.1 Substitution 


We can use substitution to convert a complicated integral into a simple one. 


Example 4.14. Consider the indefinite integral with integrand (27 + 3)1°°. We make the 
substitution 
u=2e+3 => OM i Le. es as 
da 2 
So we calculate the integral as follows: 


1 
festa dz = peo 5 du 


1 
= 5 f du 


1 


a Qv 4 101 ; 
500 | x+3) C 
We can check the result by performing the following differentiation: 
d {1 101 
Qn 4 101 = Qe 4 100 , 9=(2 100 
ag 509 (22 3) c| 509 (22 3) (29 +3)"; 


which is correct. 


Example 4.15. Suppose we have the integrand «(a + 1), Let us try the substitution 
u=x+1, sor=u-1l => —=1 ie. dr=du. 


So we have 


feet 1) ae = fo 1)u°? du 


o2 ol 
i \ 52 1 | 51, 
= 59 (2 + 1) 5 - 1) +C 
Check: 
d 1 52 1 51 51 50 50 
a eget TD pF Fe] =e ED — (et = tl)", 


which is correct. 
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Example 4.16. Consider the integrand 1/(#lnz). Let us try 
d 1 
u=Ing = “= - ie. de=cdu. 
dz 
Thus, we calculate the integral as 


1 1 
/ dx = / -xodu 
alna ru 


In |u| + C 
= In|In2z|+C. 


Check: d 
= (In| Inaz]) = 


1 
Inge a 
which is correct. 
Example 4.17. Consider the integrand 1/(1+ \/z). Let us try the substitution 
w=l+JVc => VJVe=u-l. 


Differentiating we have 


d 1 
—. = 5st? => dx = 2x2du =2u— 1)du. 
Therefore, we calculate the integral as 


laws = [2 1) du 
= 2 f (1-2) du 
= 2 f au-2 f du 


= 2u-—2In|ul| +C 
= 2(1+Vz)-2mn|1+V2z|/+C. 


Check: 
d = 1 iL 1+VJr-1 1 
2(14 2In|1 4 tC] = = = , 
dx [20 + ve) ee ve VJael+vVz) Ve+ve) 1l4+va 
Example 4.18. Consider the integrand sin(3z +1). Let us try the substitution 
w=304+1 = os Le. en 
da. 3 
So integrate as follows: 
1 
[suse +1ldx = 5 f sinudu 
: +C 
= —=cosu 
3 


1 
= ag cos(3a + 1)4+C. 


Check: 


— [5 cost + 1)4 c| = 53: sin(Bx +1) = sin(3a + 1). 
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Example 4.19. Consider the integrand sin (4) /x?. Let us try the substitution 


1 = du 1 
et — = 
xL dx x2 


So we calculate the integral as follows: 
. HE ‘é 
sin (= sin u 
/ a ve = / 2 - (—2) du 


= cosut+C 


Check: 


which is correct. 


How do we find a suitable substitution? Usually by observation and using our knowledge 
of differentiation. Or we put 


u=f(x), then “ =) = easy aide 
If we have F 
4 ipa@l= Fe), 


then we can write 


For instance, in the case 


we know that 


So we write the integral as 


1 1 
‘| : ax = f -ad(Inz) = f Satna), 
elnz alnaz Ina 


Therefore we know u = Inz will work! 


Similarly we know 


dx eg 


so the integral from example 4.19 can be written as 


[=> dx = - f sin (=) (=) u= . 
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Also, from example 4.14, we have 
1 
fe 3) Ge = 5 fe + 3)1© d(2x + 3), 


since 
d(2x + 3) 


dx 


which is the same as what we tried earlier.? 


1 
=2 dx = 5d(2x + 3) = Iletu=2r¢+3, 


°Fnd Lecture 17. 
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4.3.2 Trigonometric substitution 


Example 4.20. We know that 


1 4 
i ie = sin r+ C; 
(see pg. 38), since 
d ees | = 1 
Gy sin x) = Vilage 
Actually, we can work out this integral by a substitution like x = sinu because we know 
that 


1 — sin? u = cos? u, 


and 


dx 
du 
Thus, we calculate the integral as 


=cosu or dx=cosudu. 


cos U 


1 
dz = 
ere 1 —sin? u 


du= f du=u+C=sin 2 +C 


since 
sin~' 2 =sin7!(sinu) = u. 


Example 4.21. 


1 1 
[pert tete, since qa ttam |) aes IT 
Let us try the following 
d it 
x =tand, = cqhp = it tan’ —= dx=(1+tan?6)d0. 


So we calculate the integral as 


it 1 
—_ dz = | ———__(1+ tan? 6) dé do as 
im x / mam + tan* 0) / é+C=tan ~2+C 


Since we have 


a=tand = > tan ‘a¢=tan '(tand) = 8. 


Example 4.22. Consider the integrand 1/(1 + 227). This is similar to the previous 
example. If we try 


1 


V2r=tand => dr=—=(1+tan?6)d0, @=tan '(V2z). 


So we calculate the integral as: 


ae a SS Se EE 
1+ 22? 14+ (V2x)? 
1 1 
= / rie + tan? 0) d0 


J/2/ 1+tan? 
il 
V2 
1 
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Ceck: 
1 1 1 


d 
V21+4+ (v2z2)? es 1 + 22?’ 


= < tan7}(V/2z) 4 c| = 


which is correct. 


Example 4.23. Consider the integrand 2/(1 + 2a”), which is a variation of the above 
example. So let us try /2x = tan@. Then, we can write 


x Js tan d 1 5 
dx = 14 d 
\op - eer 75! vende 


= 5 | tnoa6 


1 f sind 
= = dé 
2 i cos @ 


1 1 
= —-~ | ——d 0). 
2 / cos 0 ions?) 
We achieve the last line from knowing that 
d 0 
(6080) _ _ sina. 
Now we use another substitution to tackle the integral we have so far, that is put u = cos 0, 
so 
Me 1 1 
aon ARON ak 
= -—-In|jul|+C 
= —=In|cos6|+C 


1 
= 5 In(1 + 2x?) + C. 


We gained the final result by using the fact that | cos6| = Vcos? 0, cos? @ = 1/(1 + tan? 0) 
and two log rules. 


This result can be gained in a quicker way, that is 


He 1 1 
——_.dr = d(1+ 2x? 
| ame kes (Ete) 


ae 


— In|ul + 


1 
7 In(1 + 22) + C, 


where we used the substitution u = 1 + 22”. 
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4.3.3 Partial fractions 


Example 4.24. Consider the integrand 1/(a?—1). We know that x? —1 = (#+1)(x—1). 


So we may write 
1 A B 


e—l 2-1 2+’ 
where A and B are numbers to be found. Re-writing the RHS of the above, we have 


Ad = SBis AAS Depa BA 
e—-1' «+1 (7-l1)(@+1) — a?-1 
Therefore, we require 
A+B = 0, 
A-B = 1, 
which has solution A = ; and B= —}. Therefore the integral becomes 


1 1 1 1 1 
> dr = d d 
/sae sfea# sfewr® 
1 1 1 1 
= d 1 d(x +1 
sf ie ) sf e ) 


1 
= =-lIn|jx-]| 
2 


= In 
2 


x—-1 
xc+l 


Example 4.25. Consider the integral 


/ v?+6e+1 
3a? + 5a — 2 


We first factories the denominator so that 327 + 52 — 2 = (3x —1)(2 +2). Also, it would 
be easier if the polynomial in the numerator was of one degree less than the denominator. 
So we manipulate the integrand as follows: 


w+6c+1 —— 3(3a?+ 18243) 
3x? + 5a — 2 3272 + 5a — 2 
_ -3(82? + 5a — 24132 +5) 
322 + 5a — 2 
1 13% +5 


3 [0° 302+52—2]” 
Now we proceed by changing the fraction in the square bracket using partial fractions, 


that is 
Be+5 A | B _ (A+3B)x+2A-B 
322 +5¢—-2 32-1 «+2 (32 — 1)(x + 2) 


Matching the coefficients on the numerator we must have 


A+3B=13 
2A-B=5 
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Solving simultaeneously we have the solution A = 4 and B = 3. So the integral becomes 


fotette if 13245 
dy = | xv 
3a? + 5a — 2 2 3a? + 5a — 2 


1 1 4 if 3 
= d | d 
al of Teste aeesta 
1 4 1 Hn ae | 
= d d 1) 4 12 
al 7 eS eee) sf ape ) 


4 
g nse 1] +In|2+2/+C. 


| 
8 


4.3.4 Integration by parts 


This is equivalent to the product rule for integration. Suppose we have two function u(x) 
and v(x). Then the product rule states 


Rearranging the above gives 


Integrating both sides we get 


‘i; uv’ dx 


[plunder — futvae 
x 
Uv - fate dx. 


So we write the rule for integration by parts as: 


dv du 
pes dx = uv — i Ue dx. (4.6) 


Given an integral whose integrand is the product of two functions, we choose one to be 


u and the other to be v’, from which we can calculate u’ and v. Plugging into the above 


equation hopefully leads to an easier integration.* 


“End Lecture 18. 
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Example 4.26. Suppose we have the integrand xe”. Let us choose 


u=a2, vi =e* u =1, oe 


Therefore, we can calculate the integral as follows: 


[ee dx = pu! dx 
uv — ue da 
= «xe* - fe dx 


= xre—e +C. 


Check: 


— [xe* — e” + C] =e” + xe” — e* = ze”. 


dx 


Example 4.27. Suppose we want to integrate Ina = 1-Inz. We choose 


/ 
u=Inzg, v=1 u , U=a. 


So we calculate the integral as 


rede = princas 
= [wae 
w~ [weds 
cine ~ f > -2dz 
x 


= «lnx—-—2+C. 


Check: F ; 
—|ehnae—2+C)]=me+2-—-l=Ing. 
x x 
Example 4.28. Suppose we want to integrate e* cos x. First let us choose 


u=cosz, vu =e* u’=-—sing, v=e’. 


So we write our integral as 


[eooserate = [wae 
= w— [weds 


= Coose+ f e*sina dr, 


Now, we have an integral similar to what we started with, so let us integrate this by parts 
too, choosing 


ui=sinz, wv! =e” a’ =cosz, t=e”. 
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So our original integral becomes 


[eeoseds — coos + f uv! de 


= Coose + wo f alvds 


= Coos +ePsina ~ fe cosrdz, 


Note, now on the RHS we have the same integral we started with. Rearranging this, we 
can make the integral the subject, i.e. 


[ereosedr = ct cosa tetsina— f €* coscedr, 


2 [6 cos. de =e" cosx +e" sing. 
So finally, we can write 
1 
fe cos xz dx = 5 le” (cosa + sinx)} + C, 


remembering the constant of integration! Check: 


d 1 
e* (cosx + sin x)] 4 c| = 5 Le" (cos x + sinx) + e*(—sinx + cosx)} = e* cosa, 


Gla | 


which is correct. 


4.4 Definite integrals 


Remark: all the techniques acquired can be applied to definite integrals. 


Example 4.29. Consider the integral 


1 
p=} (2 +1)3 dz. 
0 
Let us choose the following substitution: 
u=“2+1l = dzr=du, 
then the limits of the integral become 
et=0->u=1 and «¢=1>u=2. 


Therefore, we calculate the integral as 


1 
ay udu = —u' 
T 4 


2 


ge ono 


This is the same as finding the indefinite integral first, 
1 
ic 1)? dx = git lite, 


then imposing the limits, so 


2 
fey uw du = i@+nt+¢] 
1 


- [ile 1) +¢| 


CHAPTER 4. INTEGRATION 81 


Example 4.30. Consider the integral 


x 


2 
| cos’ x dz. 
0 


Let us choose 
‘i ° 7 
u = cos’ x, v=cosr => w=-—2cosxsinz, v=sing. 


Then the integral is calculated as 


3 
T= | cos? x cos x dx 
0 


TT 


xan 2 
= cos’ #sina ld +2 cos x sin? x dx 
0 


us 


2 
2 | cos x sin? x dx. 
0 


Now there are two ways to finish the integration. 


1. 
2 2 
PS a cos x sin* x dx 
0 
2 2 
a cos x(1 — cos“ a) dx 
0 
z a 4 
= 2 cosnde ~2 f cos” x dx 
0 0 
= 2sinz|? —2I 
= 2-21]. 
Finally we can write 
2 
P=2-2 => T= 3° 
2. 


us 


2 
I = 2 | cos x sin? x dx 
0 


us 


2 49 7 
2 sin“ xd(sin x). 
0 
Using the substitution u = sinz, then x =0 > u=0 andx= 45 +>u=l1. So, we 
have ; : 
2 2 
r=2/ uw du = <u? 
0 3 


0 3 
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4.5 Numerical integration 


Consider evaluating the definite integral 


[ f(a) dx. 


In practice, we may only know f(x) at some discrete points, and even if we know f(z), 
its antiderivative may not be expressed in terms of the functions we know, for example 


[viteae or peta. 


Since most integrals can not be done analytically, we do them numerically. We do this 
using a “geometric idea”. 


4.5.1 Trapezium method 


We want to estimate the integral of f(a) on the interval [a,b], which represents the area 
under the curve y = f(x) from a to b. 


H H H i H i i H i re 
Qa=% 1 % % Ta LH Le XW LE tg =b 


Figure 4.9: Forming trapeziums with height of the sides dictated by the curve 
y = f(x) over the interval [a,b]. 


We choose n number of pieces. Divide the interval a < x < 6 into n (equal) pieces with 
points 
a=% <2 <%2< +++ <an_-1 < an =D. 


On each piece of the interval, we build a trapezium by joining points on the curve by a 
straight line. We calculate the total area by summing all the area of the trapezia. This is 
our estimate of the integral. 


To start with, let h be the width of one piece of the interval, i.e. 
b-a 


y] 


h= 


n 
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then we have 
Lp=Xo+kh, k=0,1,2,...,n. to =a, Ln=b. 


Let us consider the trapezium based on the piece [x,_1, 2%], whose width is h. The height 
of the sides of the trapezium are f(#,_1) and f(x,). So the area is 


f(®r-1) + f(@k) 


h : 
2 
f(te-1) h f (xr) 
——— 
Le-1 Tk 


Figure 4.10: Trapezium constructed over each piece of the interval, where each piece 
has width h. 


Then the total area under the curve over [a, b] is the sum: 


pee: prod Tan} : pL + Hee} ted p Leen) Fen) 
= rao) + Fler) + (Fler) + Flea) + + lent) + Fn) 
= A rao) + 2(Fler) + Flos) +--+ Flen1)) + Fen) 
Aa 
=} flo0) +297 flax) + Han) 


We can think of the sum as follows, we have the two outer sides of the first and last 
trapezium, then every trapezium in-between shares its sides with its neighbour, therefore 
we require two lots of the interior sides.° 


5nd Lecture 19. 
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Example 4.31. Using the trapezium method, estimate 


1 
1 
fe 
0 1l+2 


We choose n = 4, then 


1—0 1 
h= = =kh, k=0,1,2,3,4. 
4 4? Xk , 0, ) 3, 
Also, note that 
1 1 
f(z) T+ 24’ ie. f(x) 1+a1 


Therefore, we have 


1 
i. dx & f(r) + 2(f (21) + f(@2) + f(%3)) + f(xa)] 


1 
8 

7 tua (BE 14 2) | ] 
8 M7. in B8r 2 
0 


a | 


Figure 4.11: Numerically integrating under y = 1/(1 + x*). Dividing interval into 4 
pieces of width h = 1/4. 


This is an over-estimate of the integral since y = f(x) is convex (i.e. it curves up like a 
cup). If it were concave (i.e. curved down like a cap), then you would have an under- 
estimate. 


Example 4.32. Estimate the following integral 


ah 1 
————-dz, ie. ¢) = ———... 
i 1— 2x4 f(@) V1— x4 


However, notice that we can’t calculate f(1). This is because y = f(x) has a vertical 


asymptote at x = 1. 
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Ry 


Figure 4.12: Graph of y = 1/(V1—2*), with an asymptote at x = 1. Estimating 
integral on interval [0, 1]. 


The problem here is that at x > 1, 1 — x* — 0, rather like 1 — x. 


Since 1-24 = (1-—2)(1+2+ 27+ 2), ie. 1— 2+ contains a factor 1 — x, which makes 
f(x) become singular at 2 = 1. So we may try to get rid of it by a substitution. 


1 1 
1 1 1 
a te= | es 
i 1-4 0 Vl-aVl+aete? +23 


Let us try u = /1—a, then u? = 1-2 or g = 1— vu? and dx = 2— udu. Now, at 
z=0>u=landx=1—-u=0. Thus, the integral becomes 


ma ra —2u 
a du 
0 Vl—a4 1 U/14+1—w+(1— v2)? 4 (1 u?)3 


0 
= 2 
i, V4 — 6u? + 4u4 — ub 
: 1 
0 V4—6u2 + 4u4 — ub 


where we have applied the rule 


[soe=- free, 


i.e. if you switch the limits, the integral changes sign. 


Let us put 
1 


rs V4 — 6u2 + 4u4 — 8 


g(u) 


Choosing n = 4, we have h = i, so 


Up =kh, k=0,1,2,3,4. 
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Then, we estimate 


1 1 1 
| ——— dz = 2 | g(u) du 
o Vliet 0 


h 1 1 3 
~ 2: + g(1 
2-5 [o0+2(9(g) +9(5) +9(2)) +00] 
= 1.32. 
The exact result is 1.311..., so we have a close estimate given we only chose 4 divisions 


of the interval. 


4.6 Application of the definite integral 


4.6.1 Area bounded by curves 


As we have discussed, integrating allows us to find the area bounded by the z-axis and a 
curve y = f(z). We can extend this to find the area between different curves on the same 
axis. 


If f(z) is a non-negative function on a < x < b, then i f(x) dz is the area between the 
curves y = f(x) and y = 0 (ie. the x-axis) from a to b. That is, the region is bounded by 


y=f(r), y=0, wr=a, r=b. 


In general, the area between the curves y = f(x) and z-axis from a to b is 


| "| #@)| de. 


(a) Use y= f(a). (b) Use y = |f(«)|. (c) Use y = |a°|. 


Figure 4.13: Examples of area bounded between curves y = f(x) and y = 0, on the 
interval [a, b]. 
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Recall: 


but the shaded area is given by 


1 0 1 1 
A= | je] dx = | (-2')de+ f a? dx = =. 
eal =a 0 2 


Suppose the region is bounded above and below by two curves y = fi(x) (top) and 
y = fo(x) (bottom) from a to b, then the area of the region is 


b 
J \i@) - plo) ax, (4.7) 


y = fila) 


y = fala) 


Figure 4.14: Area bounded between two curves on the interval [a, 6]. 


Example 4.33. Find the area of the region between y = x +1 and y= 7— 2 from x = 2 
tow =5. 


BY 


Figure 4.15: Shaded area between y= x+1 and y =7-—--« over the interval [2,5]. 
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First, we need to find the point where the two curves intersect. The point should satisfy 
both equations, therefore we solve 


r=3 
y=4 


ys=o+l 


er = j@£+1=7-c => { 


We need to know the point of intersection because from the graph, it is easy to see that 
on the left of the point of intersection, y = 7 — x is above y = x + 1, whilst on the right 
y =x+1 is above y = 7— 2. So we have to be careful when employing the formula (4.7). 


Therefore, we can finally calculate the area as follows: 


A 


3 5 
(aC 2) (e+) de+ f oe) = saide 


I| 
so 
oo 

— 
aD 

| 
1) 
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Q 
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Exercise 4.1. Find the area bounded by the curves y = 2? and y = 1—2?, on the interval 
[—1, 1]. Hint: is there more than one point of intersection?® 


®Fnd Lecture 20. 
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4.6.2 Finding a distance by the integral of velocity 


If you know the velocity u(t), then the distance s as a function of time, i.e. s = s(t), is 


s(t) = fro dt, (since s‘(t) = v(t)). 


Similarly, if you know the acceleration a(t), then the velocity can be found by the integral 
of a(t), 


v(t) = paw dt (since v’(t) = a(t)) 


Example 4.34. A ball is thrown down from a tall building with an initial velocity of 
100ft/sec. Then its velocity after t seconds is given by v(t) = 32t+ 100. How far does the 
ball fall between 1 and 3 seconds of elapsed time? 


First let us write the distance as the integral of the velocity, that is 
s(t) = ico = i 32t + 100 dt = 16¢? + 100¢ + C. 


Then the distance fallen, § say, is given by 


5 = s(3) — s(1) = (16t? + 100¢ + C)|,_, — (16t? + 100¢ + C)|,_, = 328 ft. 


Notice that 


(16¢? + 100¢ + C)|,_, — (16t? + 100+ C)|,_, = (161? + 100¢ + C)|; 


3 
/ (32t + 100) dt 
1 


= e u(t) dt. 


In general, if v(t) is a velocity function, then the change in distance between t, and tg is 


s(ts) (4) = f "Oa, 


1 


and the total distance is 


NOTE: here we call 7 a dummy variable. 


Chapter 5 


Differential Equations 


Differential equations are equations involving a certain unknown function and its derivative 
(equations for functions). For example 


d 
y =x" or i 
dx 


where y is a function of x, i.e. y = y(x). We know that in this particular case, 


1 
y(x) = =2° +C, 
3 
where C' is an arbitrary constant. 


Recall exponential growth and decay: we define the relative rate of growth as y'/y. If we 
assume that the rate is constant, say A, then we have 


y=, (y= y(t). (5.1) 
This has solution y = Ae*“, since 

y(t) = AAe™* = dy(1), 
where A is an unknown constant. 


The solution Ae is called the general solution to the differential equation in (5.1), im- 
plying that we have found all possible solutions. 


To determine the constant, we need an extra piece of information. For instance, if we 
know y at the initial time, say y(0) = yo (given), then 


y(0) = Ac’ = A= yo, 
so y(t) = yoe™. 


The problem y’ = Ay, y(0) = yo is call an initial value problem (IVP), which has a unique 
solution y = yoe**. 


The order of a differential equation is the order of the highest derivative appearing in the 
equation. 


90 
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Example 5.1. 


y' + 2xy = e® - first order, 
y” + 3y’ + 4y = 0 - second order, 
(y')? + 2Iny + 4e” = 2° - first order, 


y) + 3y"” + 2y! + 4y = 10 - third order. 


5.1 First order differential equations 


91 


Here we will consider different techniques to solve first order ordinary differential equations 


(1st order ODEs). 


5.1.1 Separation of variables 


Let us consider the simplest case of a first order differential equation, that is 


dy 


dx = Fe y), 
for example 
y! + Qary = e* y’' = —-2ay+e* = f(z, y), 
(y’)? + 2Iny + 4e* = 23 y! = +4/x3 — 2iny + 4e* = f(z, y), 
yy —x=0 y' =2/y? = f(z,y). 


Suppose that f(x,y) is separable, i.e. it can be written as 


for example 


x 1 1 
LY) = =e z)=2, hly)=-—. 
f(x,y) ry rr g(x) (y) 7 
Then we have 1 
—— dy = g(x) dz, 
a) g(x) 


(5.2) 


where the LHS only depends on y and the RHS only depends on x. Integrating both sides 


we have 


| igitr= [oan 


Therefore, if we can work out the integrals, then we can obtain the general solution to 


the equation. Also, since we expect a constant of integration, we can merge the constants 
from both sides into one constant, say C, since we are integrating the equation once only.! 


1End Lecture 21. 
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Example 5.2. Consider the differential equation 
y = ry. 


We already know the solution to this equation. Now let us see how to derive it using 
separation of variables. 


taking all things relating to y to the left, and for t to the right, we have 


1 
dy = ddt. 
y 


[oda frat 
y 


hence, using what we have learnt in previous chapters we get 


Integrating both sides we have 


Iny=At+C. 
Finally, re-arranging for y, we have 
peer! a Ae. ASE, 


Example 5.3. Consider the equation 


dy 
—— SS rb 

ae Y) 

following the procedure as in the previous example, we have 


1 
—dy = xdx. 
y 


[ow= fees, 
y 


1 
= Iny=50°+C. 


Integrating both sides we have 


Taking exponentials of both sides in order to re-arrange for y, we get 


L 


2 Lo 
gsr! = Ae, Asie’ 


=e. 


We can check if this satisfies the original equation: 


dy d 12 Te 
Be A, e\eA w+. oF = ry. 
dx an (40? — 2 ae 


Example 5.4. Consider the differential equation 
yy =. 


We first write it in the form y’ = f(z, y), ie. 
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now we realises that we can apply separation of variable, so 


y'dy = «dz, 

=> [vay= | ae 
13 lo 

ce ee 


1 
3 3 
= ye ee) 


where C” is some constant (different to C, since we multiplied through by 3). Again, we 
check the solution satisfies the equation 


1 
dy _— d Bg Sp \e 
dx dx le c’) 


=, 
ue 
Example 5.5. Consider the following initial-value problem: 
dy 2 2 
i 1 0) =1. 
qe Ye"), yO) 


First, we find the general solution, note, we can use separation of variables in this example, 
so i 
i (14+ 27) dx 


1 1 
=z+-2°+C 
y 3 
—> => 1 
0 ot gue +C 
Now check that the general solution satisfies the original differential equation: 
d d 1 1+2? 
=] 3 = 3 5 = (1+2")y’. 
zo ad&\ x£+32°+C (2+ 323+C) 


Now it remains to find the constant C’, by applying the condition y(0) = 1, i.e. we put 
x=0. 


1 1 
1 


C=-1. 
0+5:-0+C C 


y(0) = 
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So the solution to the initial value problem is 


_ 1 
ae ae ee aa 
Example 5.6. Consider the initial value problem 


aly SSa75 (0) = 2: 


First, find the general solution, 


ely! = 32” 
dy 207 
=> —=32re4 
dx 
= dy =3a7 dr 
=> pera = [30% a0 
= @an 40 
= y=lIn(2?+C). 


Check: 
3x2 a 


= 32? . 
pe eae 


dy d 7 
Aa (In(a*® + C)) 


dx dz 
Recall e!™@ = a, using this, we can write 


dy = 3226" ac) = 3r2¢e7 N@*+C) — 3926-9, 


Now we apply the initial condition, 


y(0) =In(C) = 2 C =e’, 


so we have the final solution 


5.1.2 Linear first-order differential equations 


An n-th order differential equation is linear if it can be written in the form: 


y) + an—r(a)y) + dn—a(a)y +--+ +ar(e)y' + a0(x)y = f(x), (5.4) 


or 
d”y d’-ly ay dy 
aan + On—1(2) Fay + On—2(2) 7 m5 fees ay(2)— tag(x)y = f(x), (5.5) 
where a; (i = 0,1,2,...,n) and f(x) are known functions of z. 


NOTE: here we have leading order coefficient of 1. 


NOTATION: the n-th derivative of y with respect to x is written as 


(ny — ay 


y da” 
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Example 5.7. 


y' + 2y =e", 
is a first-order linear differential equation with constant coefficients, i.e. 


rol, Ge) =2- Fle) =e 


y" + 8y' + 16y = 0, 
is a second-order linear differential equation with constant coefficients, i.e. 


n=2, a(r)=8, ao(z)=16, f(x) =0. 


y" +yy' =s, 


is a second-order non-linear differential equation. 


ety’ +x2y=cosxr => y' +2x7e *y=e *cosz, 


is a first-order linear differential equation with 


n=1, ao(x)=2e", f(x) =e * cosa. 
5. d 
— = 32 + Be%, 


is a first-order linear differential equation with 
Ro 1; cagt) =3,.. FHSre* 


Definition 5.1. If f(x) = 0 (zero function), then the linear equation (5.5) is said to be 
homogeneous; otherwise, we say the equation is non-homogenous. 


Example 5.8. 


y’ +2y =e” - non-homogeneous, 


yo +2y=0  - homogeneous. 


A first-order linear equation has the form 


oH + g(a)y = ple). (5.6) 


If p(x) = 0, then we can try to solve the differential equation by separation of variables. 


In general, we solve as follows. Let 


then 
el (2) = et q(x) dat 


is called the integrating factor. We will use this integrating factor to derive the general 
solution to the first-order linear differential equation.” 


?End Lecture 22. 
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Consider the general case for a first-order differential equation given by (5.6). First let us 
multiply both sides of the equation by e2), 


d 
cE + FP g(a)y = eT p(n). (5.7) 


Now let us consider the derivative of e7™y, 


di, ta), — .T(e)4¥ ,. 4) 72) 
AG y) = e aie yz le ) 
d d 
— ,T(«) T(x) © ¢p 
TOD + yeT © (T(0)) 
zx d zx 
ar a ye" ™q(a) 
= eT )p(z) 


Here we have put 7’(x) = q(x) and applied (5.7). So we have 


d x x 
(ey) = ep (a), 


Therefore, integrating both sides we have 


Sh) gs: ul Tp a\de. 


or 


y= eT) f M(x) dx. 


we have shown that if y’ + q(x)y = p(x) and T(x) = f q(x) dz, then 


y= ee) f Me p(a) dx. (5.8) 


Note, we expect a constant when we complete the integration above. 


Example 5.9. Consider the differential equation 


note that f(x,y) =x — (y/x) can’t be separated. So we put 
1 1 
qaya 5, palace => Tx) = [oo dz = / —dzx = Ing. 
x x 
Then we have the solution 


1/1 1 C 
goes [ htcdsa eM? | tds = ee +C| = a7 +—. 
x13 3 x 


Example 5.10. Consider 


So we put 


(yen aye. SS TiS fae) d= [ou = so. 
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Then the solution is 
y=e™ pet's dx = e723” / er” d (52) | =e” jea*” 


y=1+Ce2™, 
NOTE: this example could have been done using separation of variables. 


Example 5.11. Solve the initial-value problem: 
y=yt2, y(0)=1, 


so we put 


Therefore, the solution is 


y=e" i. e* x7 da. 


We calculate the integral using integration by parts, 
ecw de = —g%e?+ / Que” dx 
—x%e-* — Qe? + / 2e * dx 


= ere [-2? — 22] — a) 


Ser (a? + 2x + 2) +C. 


Hence, the solution to the differential equation is 


y = —(a? + 2 + 2) + Ce?. 


It remains to use the initial condition to find C, ice. 


GOS I4 Clem C =3, 


so the final solution is 


In the previous three examples, we gained the following results: 


Ex. 5.9. 1 1 
yYt-y=2, y=ar?+c-.-, 
x 3 x 
Ex. 5.10. 
y tay=2, y=1+C-e 27, 
Ex. 5.11. 
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These examples have something very important in common, that is the solutions have the 
following form 

y = f(a) + Cg(2), 
with explicit functions f and g. Here y = f(x) is a particular solution (take C' = 0) of the 
non-homogeneous equation, and y = g(x) is a solution of the corresponding homogeneous 
equation. For example 


Ex. 5.9. 
1 2 1 1 
if y= ars then y/ 4 y= goto -3r =o, 
1 1 1 i oll 
if y=-, theny’+—y= got a 
x x x oe 
Ex. 5.10. 


if y=1, theny’+ay=0+2=2, 


2 (1 
if y=e 2", then y +ay=e (5-22) +2ed* =o. 


These examples reveal an intrinsic structure of the general solution of a linear differential 
equation. We solve the equation by finding the solution of its homogeneous equation, and 
a particular solution to the non-homogeneous equation. 


Now we will understand how to use this method to solve a first-order linear differential 
equation with constant coefficients: 


y +Ay=p(x), A is constant. 
We know that the general solution is 


y(x) = f(z) 5 Cg() 
AH “~~ 


particular integral (P.I.) complementary function (C.F.) 


and 
f+rAf =ple), go +Ag=0, 


where C' is the constant of integration to be found. We start by building g. So we need 


to solve 
g =Ag =9, 


solution: 
g=Ce. 


Therefore, the general to y' + Ay = p(z) is 
y= f(z) +Ce™, 


where f is a particular solution (depending on p(x)). In what follows, we shall find f(z) 
for certain kinds of function p(x).° 


2nd Lecture 23. 
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We will consider three different kinds of options for p(x): 


1. Polynomial, 
2. trigonometric functions sin or cos, 
3. exponential function: 


(a) p has the form e™ where r 4 —A, 
(b) p has the form e’™ where r = —A. 


Example 5.12. Consider the differential equation 
yt+y=, 


so we have 
A=1, p(x) =, solution: y = f(x) + Ce’, 


i.e. we see that the C.F. from the homogeneous equation is g = Ce”. Now, here f 
should be a polynomial with a degree of one, since p(x) = x. So we try the most general 
first order polynomial, f(x) = ax + 6, and so f’(x) = a. Substituting y = f(x) into the 
differential equation we have that 


fi+f=at+ar+b=ar+(a+b)=2, 


Thus, comparing coefficients from the LHS and RHS we must have that 


a=1 a=1 
eet ~x = sa ay 


so f(z) =x —1 isa solution. Therefore, the general solution to the original equation is 
y(x) =x2-1+Ce™. 
Example 5.13. Consider the differential equation 
y +2y=a", A=2, p(x) =2?. 
We can easily work out the C.F., the general solution will take the form 
y(x) = f(v)+Ce™. 


For the particular integral, we try f(x) = ax? + ba +c, since p(x) = 27. So, we have 
& 


f' (x) = 2ax + b. Substituting into the differential equation, we have 


f' +2f =2axr +b + 2ax? + 2bx + 2c = 2ax? + (2a + 2b)z + (b+ 2c) = ex”. 


2a=1 a= 

- mea _ il g 1 1 
a+b=0 => b=—-5 => f(z) = 52 ae 
b+2c=0 c=} 


and the general solution is 


CHAPTER 5. DIFFERENTIAL EQUATIONS 100 


Example 5.14. Consider the differential equation 
y +y=sin2z, A=1, p(x) =sin2z. 
The general solution will have the form 
ylo) = f(w) +Cen*. 


To find the particular integral, we try f(x) = asin2xz + bcos2z, since p(x) = sin 2z. 
Therefore f’(a) = 2acos 2x — 2bsin 2x. Substituting into the differential equation tells us 


f' +f = 2acos 2x — 2bsin 2x + asin 2x +bcos 2x = (a— 2b) sin 2x + (2a+b) cos 2x = sin 2z. 
Comparing coefficient from LHS and RHS then 


a—2b=1 


= %55 > fla) = Zsin2x — 2 cos2 
2a+b=0 pees? Di se Ro or Oa 


5 


So the general solution is 


1 2 
y(x) = =sin 2x cos2z + Ce™*. 
5 5 
Example 5.15. Consider the differential equation 
yep Se. Nel. eaee*. Kao SS 1. 
The general solution takes the form 
y(a) = f(a) +Ce™. 
2% 


For the particular integral, we shall try f(x) = ae?*, since p(x) = e?” and f’(x) = 2ae”. 


Substituting into the differential equation gives 
f af =Saer aoe" = Soe Se 
Comparing coefficients gives a = $ thus the general solution is 
y(z) = so +Ce*. 
Example 5.16. Consider the differential equation 
y+y=e*, A=1, ple) =e™. 
The general solution has the form 
ylo) = f(a) +Cer*. 
Now, if we try f(x) = ae~*, then 
f'+f =-ae~* +ae* =0Fe™%, 
so it doesn’t work since r = —1 = —X. But we may try 
f(x) =axe~*, then f’(x) = ae-* — aze™®. 


Substituting into the differential equation gives 


f'(x) + f(x) = —ave™® + ae~* + axe™*” =ae* =e", 


therefore, comparing coefficients, we have a = 1, so the general solution is 


y(z) =ae *+Ce* =e *(44+C). 
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In general, for 
Ax 


y tr»ay=e*, 

the general solution is of the form 

y(x) = f(z) + Ce. 
For the particular integral we try f(x) = are~>*, then 

f' +Af = ae — Naxe~** + Aare” = ae” =e, 
therefore, comparing coefficient, we have a = 1. So the general solution is 
y(a) = ae" + Ce* Se P(e £0). 

This can be easily generalised for the case 


y! +Ay = be~**, ~—- and b are constants. 


5.1.3 Second-order linear differential equations with constant coefficients 


y +ry +sy =p(2), (5.9) 
where r and s are constant. The general solution to (5.9) has the form 
y(a) = f(@) + g(a), 


where y = f(x) is a particular integral from the solution to (5.9) and y = g(x) is the 
complementary function from the solution to the homogeneous equation 


y" +ry' + sy =0. (5.10) 


Solving (5.10) to build y = g(x) 


Two functions (say, two solutions to (5.10)) are said to be independent if one is not a 
constant multiple of the other; otherwise, they are said to be dependent. For example 


e yi(x) = 1, yo(a) = x, here y, and yo are independent. 
e yi(x) =e, yo(x) =e”, if a #b, then y; and ye are independent since 


y= et = ebtt(a-b)x = eb-a)oa, 


ie. if a 4b, then the multiple is not a constant. 


ax+b 
b 


, a and 6 are constant, then y; and yo are dependent, 
is constant.* 


e y(t) = ce, w(x) = e 
since yo(x) = ey; (x), € 


“End Lecture 24. 
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If y1(%) and y2(x) are independent and they are the solutions to (5.10) then the general 
solution to (5.10) can be written in the form 


y(x) = Ciyi(x) + Cryo(z), (5.11) 


where C, and C2 are constants. In other words, if you have two independent solutions, 
then you can represent any other solution in terms of these two solutions. 


Now, we verify that (5.11) is the solution to (5.10): 


yl +ry tsy = Cry{ + Ciry, + Cisy, + Cayg + Cary + Cosye 
= Ci(y{ +ryy + syi) +Co (ys + rys + sya) 
Re ——— ~~ lee ——S 
=0 =0 


Therefore we need to find two independent solutions to (5.10). But how do we find them? 
Let us assume that the solutions are of the form 


Then 
y! = er. y" = ver: 


and so substituting into (5.10) we have 
y” + ry! + sy = A7e** + rre** + se*” = e** (1? +r +8) =O, 
which tells us that if A is a root of the equation 
M+rrA+5=0, (5.12) 


then y = e** will be a solution to (5.10). Equation (5.12) is called the auxiliary equation 
or characteristic equation and is quadratic, so we expect two solutions for \ and thus two 
independent solutions to (5.10). 


Example 5.17. Consider the second-order differential equation 
y” + 3y' + 2y =0. 
If we put y = e>”, we see that the auxiliary equation is 
M4 3A42=0 — > (AF1)(A42) =0. 


This has two roots: 


Therefore we have two solutions e~* and e~?*, and they are independent. So the general 
solution is 
y(x) = Cie * + Coe”. 


In general, we know that the roots of (5.12) are given by 


A= ae A =r? —As. 


There are three cases, depending on the value of A. 
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(a) A > 0: two distinct real roots (see Fig. 5.1(a)), 


—r+VA —r—VJA 
2 : 2 


AL = 


We have two independent solutions to the differential equation (5.10), 


Aix Ax 
9 


e€ and e 


so 
g(x) = Ce + Cre™*, 


is the general solution to (5.10), ie. the C.F. of (5.9). 


: Ki 
r2 M1 x 


Mt » e mn 


(a) Two roots. (b) One roots. (c) No roots (complex). 


Figure 5.1: Different options for the curve y= 7 +r +s when solving the equation 
M+rvA+s=0. 


(b) A =0, (5.12) has one root which is real (see Fig. 5.1(b)), given by 


4 . 


So e** is one solution to (5.10). We need to find another solution, which is inde- 


pendent of e*!*. So we try 


y= rer, 


then 
yl =e" + \yre™® and y” = dye™* + Aye + APae™® = N2xe™® 4+ 2 e™!*. 
Hence 


y +ry tsy = Nae* +2)e™* + re + pyre! + srer* 
(A? + 7d1 +8) ve™® + (241 +r) e%? 
=0 —0 
= 0. 


That is, re*!” is a solution to (5.10). Also, e*!” and xe*!” are independent. There- 
fore, in this case the general solution of (5.10) is 


g(x) = Cye™* + Core™”. 
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Example 5.18. Consider the differential equation 
y" + 2y' +y =0. 


The auxiliary equation is 
PEO eit: ety 0S po, 
Ai” and xe™!® are two independent solutions. So the 


where A, is the only root. But e 


general solution is 
=e *(CL + Coz). 


y(xz) = Cle* + Cone” 


(c) A <0, there is no real solution to (5.12) (see Fig. 5.1(c)). But it is still possible to 
find two independent solutions, these will be complex roots given by 


== Se See ST —1 —A 
je ea r (—1)( Pee r 5 =, A eae 
where —A > 0 and 
¢=4/ 1, a= - B= 5V-A= 5 V4 ieee 


In this case, the solutions are e(¢+#)* and e(¢—#9)*, so we have the general solution 


g(x) = Cyelaribe 4 Cy e(o-tB)e — pat (Cre fi Cxe"#*) . 
But recall Euler’s formula e’® = cos@+isin 6, which means we can write the general 


solution as 
g(x) =e ((C1 + C2) cos 6x + (C1 _ C2) sin Ba) =e (Ci cos Ga + Co sin Br), 


where e°* cos Ga and e°” sin Bx are two independent solutions. 
Let us verify 
y =e cos Bx 


is a solution of (5.10). We have 
y’ = ae cos Ba — Be sin Bx, 


= a’e™ cos Bx — abe sin Bx — aBe™ sin Bx — Be cos Bx 
(a? — 6?)e* cos Bx — 2ahe sin Bx. 


Plugging these into equation (5.10), we have 


y" +ry! + sy = (a? — B?)e°* cos Bx — 2ahe sin Br+ 
+ rae cos Bx — rBe™ sin Bx + se cos Bx, 


which can be written as 
yl! +ry' + sy = (a? —~B+ar+ s) e* cos Bx — (2a +1) Be sin Bx = 0. 
—s_ —_ 
=0 


Similarly we can show y = e®” sin 3z also satisfies (5.10), thus g(x) above is a general 


solution to (5.10). 
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Example 5.19. Consider the differential equation 
y” — 6y' + 13y = 0. 


The auxiliary equation is 
7 — 6A +13 =0, 


which has 


r=-6, s=13, A=r*?—4s=36-—52=-16<0, 


a=-7=3, B=5/-(-16) =2, 


therefore e®” cos 2x and e®” sin 2x are two independent solutions, so 


y= eC; cos 2a + C2 sin 22), 


is the general solution.® 


5Fnd Lecture 25. 
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Finding a particular solution to y” + ry’ + sy = p(x) 


It depends on the function p(«). We only consider three kinds of p(x): 


1. Polynomials, 
2. trigonometric functions, 
3. the exponential functions: 


i. e#*, ys is not a root of of the homogeneous equation y” + ry’ + sy = 0, 


ii. e4”, wis a root of the homogenous equation. 
Example 5.20. Find the general solution to the differential equation 
y! +2y +y = 7. 


Recall, the general solution takes the form y = f(x) + g(x). First we find the general 
solution to the homogeneous equation 


y" + 2y" + y =0, 
i.e. we seek the complimentary function (C.F.) y = g(x). The auxiliary equation is 


PONT | pee Mast 


i.e. it has only one root. So the C.F. is 


g = Cie" + Core™” = (Ci + Cox)e™”. 


Second, we find the particular integral (P.I.), we try 


f=azr?+br+-a, since p(x) = 2. 


so we have 
f=2ar+b, f” =2a. 


Substituting y = f(x) into the differential equation gives 


b) +ax* + br +e 


b)a + 2a + 26+ ¢ 


fU +2f' +f = 2a+2(2azr 4 
4a 4 


( 


I 
a 
8 


a=1 a=1 
4a+b=0 = b=—4 = f(x) =27-42+6, 
2a+2b+c=0 c=6 


Finally, we can write the general solution as 


y(x) =7?—4¢+64 (Ci + Cox)e~”. 
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Example 5.21. Solve the following initial-value problem: 
y"—2y’+y=sine, y(0)=—-2, y/(0)=2. 


Note, we have two conditions here because second order differential equations have two 
constants of integration to be found. The auxiliary equation for this problem is 


YAS Bo a1 So ye, 
i.e. we have a repeated root and hence, the C.F. is 
g(x) = (Cy + Coxe”. 


To find the P.I. we try 
f =asinz + bcosz, 


since p(x) = sinx. Then we have 
f' =acosx—bdsinz, f” = —asinz — bcosz. 


Substituting into the differential equation we have 


f’ —2f'+f = —asinz — bcosz —2acosx + 2bsinz + asinz + bcosx 
(—a+ 2b+ a)sing + (—b— 2a + b) cosa 


= 2bsinx —2acosx 


= sing. 


Comparing coefficients, we have 


1 
a=0, 5b 5 f 5 COS. 


Therefore the general solution to the initial-value problem is 


1 
y(x) = 5 COSz + (Cy + Cox)e*. 


In order to find the unknown constants C; and C using the initial conditions, we need to 
find y/(x), so we differentiate the above to give 


1 1 
yf (2)= 5 sine + Cye” + (Cy + Cox)e” = 5 sina + (Cy + C2 + Cox)e”. 


Put x = 0, so 


9 
Ci =—--x, C2 =2-C1=2+52= 5. 


Finally, the solution to the initial value problem is 


1 1 
y(xz) = 5 cose + ~e* (9x — 5). 
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Example 5.22. Find the general solution of the following differential equation, 
y” + 4y' + 3y = 5e*. 
The auxiliary equation is 
MW+443=0 — = (A4F1)(A43) =0. 


Hence, this has two distinct real roots, namely 4; = —1, A2 = —3. So the C.F. (from the 
homogeneous equation) is given by 


g(x) =Cie"* + Coes: 


To find the P.I. we try 
f =ae™, 


since p(a) = 5e*x. Differentiating, we have 
if = 4q* f' = 16a”. 
Substituting into the differential equation we see that 
fl" +4f' +3f = 16ae*” + 16ae** + 3ae*” = 35ae*” = 5e*. 


Therefore we have a = 1/7 and so the P.I. is 


f= qt 


Finally, the general solution is 


Example 5.23. Solve the initial-value problem given by 
y+ 4y'+3y=e", y(0)=0, y'(0) =0. 


We know from example 5.22 that the C.F. to the homogeneous equation y” + 4y' + 3y = 0 
is 
g = Cie* + Coe **”, 


ie. Ay = —1 and Ag = —3. 


Now let us find the P.L., if we tr = ae *, we know that it wouldn’t work since ae” is 
y 


actually a solution to the homogeneous equation y” + 4y’ + 3y = 0. Therefore, we try 


f=are”, 


thus 


fi’ =ae*—azve*, and f” =—2ae*+azre™*. 


Substituting into the differential equation we have 


x 


f" +4f'+3f = -—2ae~* + axe~* + 4ae~* — 4axe~* + 3axre™* 


= 2ae* 
=e”, 
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Therefore we must have a = 1/2 and so the general solution to the differential equations 
is 1 
ie) = gre” Ole + Oe 


Differentiating the general solution we have 


va \=—e* ner = Che eS 3058, 
Putting « = 0, we have (from the initial conditions) 


y(0) = C1 + C2 =0 Ci = 
y{(0)=4-Ci-3@@=0 f 


thus 


1 1 i 


iz, Bo 4 


y(x) = give AS T ig ’ 


is the solution to the initial-value problem.® 


®Fnd Lecture 26. 
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Example 5.24. Find the general solution to the equation 
y+ 2y’ +y =2e*. 
The general solution takes the form y(#) = f(x) + g(2). 
Find the C.F.: The auxiliary equation for the above differential equation is 


W4+2441=0 — (A411)? =0 A =, 


i.e. we have a repeated root so 
g(x) = Cie * + Core™”. 
Here e * and xe * are two independent solutions to the homogeneous equation 


y! +2y’ +y=0. 


Find a P.I.: We have to try 
f =az’e*, 


since e~” and xe~* can’t be the solution to the original differential equation as they satisfy 
the homogeneous equation. So we work out the derivatives 


f' = 2are~* — ax?e~*, 


f"” = 2ae~* — 2are~* — Qaxe~* + ax*e~* = axe~* — daxe~* + 2ae~*. 
Substituting y = f(x) into the differential equation, we have 


fo +2fi+f = axe -* — daxe-* + 2ae* + dare * — 2ax*e-* + ax*e* 


= 2ae* 


= 2e%, 
therefore we have a = 1. So finally, we have the general solution 


y(x) = (C1 + Con + x*)e™. 


5.2 Simple Harmonic Motion (SHM) 


SHM is essentially standard trigonometric oscillation at a single frequency, for example a 
pendulum. 


An ideal pendulum consists of a weightless rod of length | attached at one end to a 
frictionless hinge and supporting a body of mass m at the other end. We describe the 
motion in terms of angle 6, made by the rod and the vertical. 
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mg sin 0 


mg 


Figure 5.2: Sketch of a pendulum of length | with a mass m, displaying the forces 
acting on the mass resolved in the tangential direction relative to the motion. 


Using Newton’s second law of motion Ff’ = ma, we have the differential equation 
—mgsiné = mlé, 


which describes the motion of the mass m, where the RHS is the tangential acceleration 
and the LHS is the tangential component of gravitation force. 


NOTATION: 6 = d6/dt and 6 = d?6/dt?. 
We re-write the equation as 


6+wsind = 0, w= 7 
This is a nonlinear equation, and we can not solve it analytically. 


Approximation: if @ is small, then sin@ ~ 6, and in this situation we have an approximate 
equation given by 
6+wd =0. 


We solve the equation for @(t). Here we have r = 0, s = w* and A = —4w? < 0. The 
auxiliary equation is 


vM+w?=0 Maw? = A=iw, w= -iv, 


where i = /—1. So 
r 1 -— 
a 5 and £6 5M Ww, 
therefore we have e™ cos Bt = coswt and e™ sin Gt = sinwt. Hence, the general solution is 
6(t) = Acoswt + Bsinwt. 


Differentiating we have 


O(t) = —Awsinwt + Bw coswt 


6(t) = —Aw* coswt — Bw? sin wt. 
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We can verify 6(t) satisfies the original differential equation as 


6(t) = —Aw? coswt — Bu” sinwt = —w* (Acoswt + Bsinwt.) = —w?6(t). 
eee SS 
=6(t) 


The solution 6(t) can be written as 


A B 
2s, Yad 2 ; 
BS aioe (aE ae) 
V A? + B?(sin dcoswt + cos dsin wt) 
= Rsin(wt+ ¢). 


cos wt 4 


R= vA? +4 B? 


A 


Figure 5.3: Using Pythagarus’ theorem to write the constants A and B in terms of 
the phase angle ¢. 


R - amplitude of the motion. 
o - phase angle, i.e. the amount of shift. 


w = 1/g/l- the natural frequency, i.e. the number of complete oscillations per unit 
time. 


T = 27/w - period, the time taken for a complete cycle (two complete swings). T 
depends on the length of the pendulum, but doesn’t depend on the mass and initial 
conditions. 


Figure 5.4: Graph showing the change in @ over time t, displaying oscillations of 
period T. 
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If the pendulum is initially at rest, i.e. 0(0) = 0, 6(0) = 0, then 


6(0) = A=0, 6(0)= Bw =0 B=0 0(t) = 0, 


i.e. the pendulum will remain at test for all time t. 


If the pendulum is displaced by an angle 69 and released, then 0(0) = 09 and 6(0) = 0, so 


6(0) = A=%, 6(0) = Bu =0 B=0, 


therefore 
6(t) =O9coswt = > |O(t)| < 6. 


That is, if the displaced angle 09 at initial time is small, then the small angle approximation 
makes sense. 


The solution tells us the oscillation, once started, goes forever. But in reality, the frictional 
force and air resistance would eventually bring the pendulum to a rest.” 


"End Lecture 27. 
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5.3 Solving initial-value problems numerically: 
Euler’s method 


Most differential equations can not be solved analytically, so we try to solve them numer- 
ically. 


Suppose we have an initial-value problem: 


“au = f(x,y), y(a) = yo. 


We want to find the solution y(x) numerically on the interval [a, 6]. 
First we divide [a,b] into N subintervals by the points 


a=% <4 <%2< +++ << ap << aN =). 


Xo Ly x2 UN-2 LN-1 LN 
—#_+—_—1-Pd4HY+1.+4—1_1+—"+ 
a b 


If the subintervals have equal length, say h, then 


i= 
Lp =at+kh, h= (step size), k=0,1,...,N. 


Assume that y(x) is the solution we want, then 


 y(2) = flava), va) =o. 


Integrating both sides on the subinterval [x,, 7,41], we have 


[over = [fea 


k k 


Considering the LHS, we have 
LHS = y(x)|[2h** = y(we41) — y(ze), 


thus 


Uk+1 


y(@e41) — Y(@k) = / f(a, y(a)) dx. 


Lk 


Let g(x) = f(z, y(x)), then 
RHS = / - g(x) dx, 


k 


representing the area under the curve y = g(x), between x, and 2441. 
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i i 
Xk, Tk+1 


Figure 5.5: The integral which represents the area under the curve y = g(x) is 
approximated using rectangles. Error depends on the width of the rectangles, i.e. the 
number of sub-intervals. 


If we use the area of the rectangle 


(te+1 — Te) g(rE) = hg(re) = hf (xe, y(ee)), 


to approximate the area, we have 
y(@e41) — Y(@e) © hf (we, yw): 
k=0: y(a1)—yo © h(F(#0, yo) => y(a1) © yothf (ao, y1) = y1, an approximation 
to y(x1). 


k = 1: y(x2) — y(a1) © A(f(ei,y(e1)) => (2) & y(wi) + Af(e1,y(@1)) © 
yi thf (a1, y1) = yz, an approximation to y(22). 


In general, we have 
Yeti = Vk + hf (xx, ye), k=0,1,...N—1, 


where y, is an approximation to y(x,). This is a difference equation and we can solve it 
iteratively. This method is based on the above formula, and is called Euler’s method. 


Example 5.25. Estimate y(1), where y(zx) satisfies the initial-value problem: 


dy _ 


= =i 
Tn y(0) 


We know the exact solution is 
yfaj)=e", => y(1) =e & 2.71828. 
Now we apply Euler’s method to the problem. We have 


f(x,y) =. 


First, we take N = 5, then h = (1 — 0)/5 = 0.2. 
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wi) weal x2 x3 x4 x5 

+—t+—_ M1 ++ 

0 dD 2 3 4 1 
5 5 3) 


yo = y(0)=1 
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y = yothf(xo,yo) =14+0.2x 1=1.2 
yo = yr thf(r1,y:) = 1.2 +0.2 x 1.2 = (1.2)? 

yz = yothf(ae, ye) = yo + hye = yo(1 +h) = (1.2)? x 1.2 = (1.2)8 
t= LO 

ys = (1.2)? © 2.48832. 


For N = 5, we have 


error = € — ys = 2.71828 — 2.48832 = 0.22996. 


Now, we double the number of subintervals: N = 10, h = 0.1 then we need 10 steps to 


reach x19 = 1. 
yio = (1.1)1° = 2.59374, 


then we have 
error = 2.71828 — 2.59374 = 0.12454. 


For N = 20, h = 0.05 and so 


yoo = (1.05)?° = 2.65330, error = 0.0650. 


For N = 40, h = 0.025 and so 


yoo = (1.025)*? = 2.68506, error = 0.0332. 


Euler’s method is first order, i.e. the error behaves like O(h). 


In general, if h = 1/N, then 


1 
Y= yothy=(1+h)y=1t+h=1+ 


1 2 
y2 = n+ hn = (14+ An = (145) 


Thus 


Actually, 


CHAPTER 5. DIFFERENTIAL EQUATIONS 


The source of the errors when approximating the function come from 


1. discretisation error, 


2. round-off error.® 


5nd Lecture 28. End of course. 
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